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Abstract 

There has been significant progress in the understanding of the r-process over the last ten 
years. The conditions required for this process have been examined in terms of the parameters 
for adiabatic expansion from high temperature and density. There have been many developments 
regarding core-collapse supernova and neutron star merger models of the r-process. Meteoritic 
data and observations of metal-poor stars have demonstrated the diversity of r-process sources. 
Stellar observations have also found some regularity in r-process abundance patterns and large 
dispersions in r-process abundances at low metallicities. This review summarizes the recent results 
from parametric studies, astrophysical models, and observational studies of the r-process. The 
interplay between nuclear physics and astrophysics is emphasized. Some suggestions for future 
theoretical, experimental, and observational studies of the r-process are given. 
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1 Introduction 

The basic framework for understanding the origin of the elements was summarized in the classic works 
by Burbidge et al. PP and Cameron [2j in 1957 (see [3j for a more up-to-date review). The crucial 
observational basis for developing this framework was the abundance distribution of stable and long- 
lived nuclei in the solar system. The slow (s) and the rapid (r) neutron-capture processes were proposed 
as the dominant mechanisms for producing the elements heavier than Fe. The nuclear systematics of 
these two processes provides a simple and beautiful explanation for the peaks at mass numbers A = 80, 
88, 130, 138, 195, and 208, respectively, in the solar abundance distribution (see Fig. 1). The work by 
Seeger et al. [5J in 1965 gave a comprehensive discussion of the s-process and the r-process based on the 
underlying nuclear physics. More recent reviews of the s-process can be found in IHI^IHI- This review 
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Mass Number (A) 

Figure 1: The logarithm of the solar abundances (normalized so that the elemental solar 
abundance of Si is 10^) as a function of mass number A. For clarity, only the abundances of 
the nuclei with even proton and neutron numbers are shown. When there are two or more 
stable nuclei with the same mass number, the abundance of the dominant isobar is shown. 
See |4j for the full set of abundances. The peaks a.t A = 80, 88, 130, 138, 195, and 208 
discussed in the text can be clearly seen (the abundance of ^^^Pt, which defines the peak at 
A = 195, is essentially the same as the abundance shown for A = 194). It is interesting to 
note that the abundances of ^^^Sn and ^^'^Sn (with the magic proton number Z = 50) are 
comparable to those in the peak ai A = 130. 

focuses on the progress in the understanding of the r-process since the last major review by Cowan et 
al. in] in 1991. After a basic introduction, the astrophysical conditions for a successful r-process are 
discussed in Sec. 2. Recent developments regarding core-collapse supernova and neutron star merger 
models of the r-process are reviewed in Sec. 3. Data from studies of meteorites and metal-poor stars 
and their implications for the r-process are discussed in Sec. 4. Some suggestions for future theoretical, 
experimental, and observational studies of the r-process are given in Sec. 5. For those readers who are 
interested in some particular aspects of the rich physics and astrophysics associated with the r-process, 
a good place to start is Sec. 5, where a brief summary is also provided. Together with the table of 
contents, this summary may serve as a guide in finding the topics of special interest that are covered 
here. 

1.1 The s-process 

In the s-process, neutron capture is slow compared with P decay so that an unstable nucleus produced 
by neutron capture /3-decays to its stable daughter before it captures another neutron. Consequently, 
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the s-process proceeds close to the /3-stabihty hne. The number of the stable nuclei with mass number 
A produced in an s-process event, i.e., the yield Ys{A), is governed by 




cross section (appropriate for the conditions in the s-process environment) for the stable nucleus with 
mass number A. When an equilibrium is reached. 



which means that the s-process yield of a stable nucleus is inversely proportional to its neutron-capture 
cross section. As the stable nuclei ^®Sr, ^^^Ba, and ^°^Pb with the magic neutron numbers N = 50, 
82, and 126, respectively, have extremely small neutron-capture cross sections, the s-process produces 
abundance peaks at these nuclei, which are observed in the solar abundance distribution (see Fig. 1). 

The s-process encounters the magic neutron numbers N = 50, 82, and 126 at the stable nuclei with 
proton numbers Z = 38, 56, and 82, respectively. As neutron capture is rapid compared with (3 decay 
during the r-process, this process proceeds on the neutron-rich side of the /3-stability line and encounters 
the same magic neutron numbers A^" = 50, 82, and 126 at unstable neutron-rich progenitor nuclei with 
smaller proton numbers Z ~ 30, 48, and 69 corresponding to ^4 ~ 80, 130, and 195, respectively. The 
production of these progenitor nuclei is favored by the r-process due to their relative stability associated 
with the magic neutron numbers. The successive /? decay of these nuclei following the cessation of rapid 
neutron capture then gives rise to the peaks at ^^Se, ^^°Te, and ^^^Pt observed in the solar abundance 
distribution (see Fig. 1). 

It is useful to decompose the solar abundances into the separate contributions from the s-process 
and the r-process, respectively. This decomposition is trivial for the nuclei that can only be produced 
by the s-proccss or the r-process (i.e., the s-only or r-only nuclei). For example, the stable nucleus 
^^^Ba has a more neutron-rich stable isobar ^^^Xe. The successive (3 decay of the r-process progenitor 
nucleus with A — 136 stops at ^^^Xe and does not contribute to the abundance of ^^^Ba. Thus, 
^^^Ba is an s-only nucleus. On tlic other hand, the unstable nucleus ^^^Te is sandwiched by its two 
stable isotopes ^^^Te and -"^^^Tc. The s-process follows the j3 decay of ^^'^Te and bypasses ^^^Tc. Thus, 
^^®Te is an r-only nucleus. In addition, as the region immediately beyond the heaviest stable nucleus 
^°^Bi is populated by very short-lived nuclei, the long-lived nuclei beyond ^°^Bi such as ^^^Th, ^^^U, 
and ^^^U cannot be produced by the s-process and must be attributed to the r-process. Except for 
the s-only and r-only nuclei, a heavy nucleus generally receives contributions from both the s-process 
and the r-process. The decomposition of the solar abundances into the solar s-process and r-process 
abundances, respectively, can be accomplished by a phenomenological approach that is based on mostly 
nuclear physics and observational data or by a full approach that is based on nuclear physics, stellar 
physics, and Galactic chemical evolution. In both approaches, the goal is to calculate the solar s-process 
abundances directly from theory and then to obtain the solar r-process abundances by subtracting the 
s-process contributions from the total solar abundances. 

The main nuclear physical input for the s-process is the neutron-capture cross sections. The astro- 
physical input is the intensity and the duration of the neutron flux in individual s-process events, which 
can be determined from stellar physics and Galactic chemical evolution in the full approach. In the 
phenomenological approach, the astrophysical input is parameterized as a distribution of the integrated 
exposure to the neutron flux. A given neutron-exposure distribution specifies a function 



where Ns{A) is the s-process abundance of the stable nucleus with mass number A that is produced 
by this neutron-exposure distribution. The neutron-capture cross sections and the solar abundances of 
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the s-only nuclei in different mass regions can be used to select the neutron-exposure distribution that 
gives the function /q ,,(^4) as defined in Eq. Q but for the solar s-process abundance Nq^s{A). The 
solar r-process abundance NQ^ri^) can then be obtained as 

N^,r{A) = Nq{A) - N^^M) = NeiA) - (4) 

where Nq{A) is the total solar abundance of the stable nucleus with mass number A. Figure 2a shows 
the solar r-process abundances calculated from this phenomenological approach by Arlandini et al. |lUj . 

Clearly, the solar r-process abundances obtained above are subject to possibly large uncertain- 
ties when the solar s-process fraction (3q^s{A) = Nq^s{A)/Nq{A) is close to unity. For example, the 
phenomenological approach gives /^©^^(^^^La) = 0.83 while stellar-model calculations (see below) give 
/^0,s(^^^La) = 0.62 ^lOj. This difference indicates that the uncertainty in the solar r-process abundance 
of ^^^La may be as large as a factor of 2.2. An important issue is then the accuracy of the solar s- 
process abundances as calculated from the phenomenological approach. The systematic error of this 
approach can only be assessed by the full approach that calculates and incorporates the s-process yields 
in different astrophysical environments over the Galactic history prior to the formation of the solar 
system. The phenomenological approach has identified a weak s-process, which produces the nuclei up 
to ^^Sr, and a main s-process, which produces the nuclei *®Sr and above. The weak s-process occurs in 
massive stars while the main s-process occurs in low-mass asymptotic giant branch (AGB) stars (e.g., 
[U]). Arlandini et al. ^U] calculated the s-process yields of two AGB stars with a metallicity of half 
the solar value but with masses of 1.5 and 3 Mq {Mq being the mass of the sun), respectively. 
It is encouraging that the main component of the solar s-process abundances as calculated from the 
phenomenological approach can be essentially reproduced by averaging the s-process yields of these two 
stars (see Fig. 2b). However, this should be viewed only as a significant step towards checking and 
correcting the phenomenological results by the full approach. In using the solar r-process abundances 
obtained by subtracting the s-process contributions as calculated from the phenomenological approach 
or a few sample stellar models, it is important to recognize that these results are generally reliable when 
Pq^s{A) is small, but may have substantial uncertainties when ^Q^siA) is close to unity. In addition, the 
solar r-process abundances at A = 206-209 shown in Fig. 2 may have to be significantly revised due to 
important contributions from low-metallicity AGB stars |TT] . 

Another potential problem with the solar r-process abundances obtained by subtracting the s- 
process contributions is that some nuclei may be produced by processes other than the s-process and 
the r-process. The small contributions from the so-called "p-process" were considered in [U]. However, 
it is rather difficult to estimate the contributions to the solar abundances of the nuclei above the Fe 
group but with A<90 that may come from the a-process to be discussed in Sec. 12.31 In view of this 
difficulty, the solar "r-process" abundances at A < 88 with a peak at A = 80 (possibly corresponding 
to the progenitor nuclei with = 50) will not be discussed here. 

1.2 The r-process 

The pattern of the solar r-process abundances (the solar r-pattern) obtained in connection with the 
s-process studies has played an essential role in the understanding of the r-process. As mentioned in 
Sec. 11.11 the peaks a.t A = 130 and 195 in this pattern can be accounted for by the nuclear properties 
associated with the magic neutron numbers = 82 and 126, respectively. It is convenient to consider 
the nuclear physical aspects of the r-process separately from the astrophysical aspects that treat the 
conditions in the r-process environment. Provided that an r-process occurs, the resulting yield pattern 
is mostly determined by the nuclear systematics of the interplay between neutron capture, photo- 
disintegration, f] decay, and possibly fission. For the nuclei with proton numbers Z < 80, fission may 
be ignored. Then the number of the nuclei {Z, A) with proton number Z and mass number A produced 
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Figure 2: The logarithm of the solar r-process abundances (normalized so that the elemental 
solar abundance of Si is 10^) as a function of mass number A. These abundances are 
obtained by subtracting the s-process contributions calculated from (a) the phenomenological 
approach and (b) models of two AGB stars. See (TU] for details. 
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in an r-process event, i.e., the yield Y{Z, A), is governed by 

Y{Z, A) = A - 1))Y{Z, A-l) + A^,„(Z, A + 1)Y{Z, A + l) + \p^{Z - 1, A)Y{Z - 1, A) 

+Xf3i{Z - 1, A + l)Y{Z -l,A + l) + \p2{Z -l,A + 2)Y{Z - 1, A + 2) 
+Xp^{Z - 1, A + 3)r(Z - 1, A + 3) - A))Y{Z, A) - A^,„(Z, A) 

-[A^o(^, A) + A^i(Z, A) + A^2(^, v4) + A^3(^, A), (5) 

where n„ is the neutron number density, ?T,„(f ^^^^(Z, A)) is the thermally averaged neutron-capture 
rate, A^^„(Z, A) is the photo-disintegration rate, and X^q^Z^A), \pi{Z,A), A/32(-Z^, and X^^^Z^A) are 
the rates for j3 decay followed by emission of 0, 1, 2, and 3 neutrons, respectively. Equation (0) and the 
like form an r-process reaction network. Some approximations that are either helpful in understanding 
the results from full network calculations or used in place of such calculations are discussed below. 



1.2.1 (n, 7) ^ ilin) equilibrium 



When both neutron capture and photo-disintegration occur much faster than (3 decay, a statistical 
(n, 7) # (7,''^) equilibrium is achieved and the yields Y{Z,A) and Y{Z,A + 1) of two neighboring 
isotopes satisfy 



Y{Z,A + 1) nn{vanr,iZ,A)) 



Y{Z,A) 



A^,„(Z,A + 1) 



m^kT ) 



A + iy/' G{Z,A + 1) 
2G{Z,A) 



exp 



SJZ,A + 1) 



kT 



(6) 



where T is the temperature, h is the Planck constant, k is the Boltzmann constant, m„ is the atomic 
mass unit, G{Z, A) is the nuclear partition function, and Sn{Z, A + 1) is the neutron separation energy. 
The most abundant isotope in the isotopic chain of Z has a neutron separation energy of 



niukT \ 




\W K 



3/2' 



2.79 + 0.198 




MeV, 



(7) 



which can be seen from Eq. (jH)) by setting Y{Z,A + 1) ^ Y{Z,A) and neglecting the relatively small 
differences in the nuclear partition function and the mass number. As only depends on n„ and T, 
the most abundant isotopes in different isotopic chains have approximately the same neutron separation 
energy, which is ~ 2-3 MeV for typical conditions during the r-process. Due to the odd-even effect 
caused by pairing, the most abundant isotope always has an even A^. For this reason, it may be more 
appropriate to characterize the most abundant isotope in the isotopic chain of Z by a two-neutron 
separation energy ^2n(^, A + 2) = Sn{Z, A + 2) + A^V)^2Sl^. 

The conditions for (n, 7) ^ (7, n) equilibrium were first examined in ^31 based on steady- flow 
calculations (see Sec. I1.2.2|) and were found to be T > 2 x 10^ K and n„ > 10^° cm~^. Later studies 
emphasized the effects of the nuclear physical input, such as the masses of neutron-rich nuclei 
far from stability, on the conditions for (n, 7) ^ (7, n) equilibrium. Essentially all of the nuclear 
physical input for the r-process is calculated from theory (see |9j for a review) and different calculations 
sometimes give very different results (e.g., f2lElE])- The uncertainties in the theoretical results can 
be reduced by making use of the existing but limited experimental knowledge on neutron-rich nuclei far 
from stability (e.g., JHl)- Based on full network calculations with two different sets of nuclear physical 
input, it was found that (n, 7) ^ (7, n) equilibrium is obtained for n„ > 10^° cm~^ at T = 2 x 10^ K and 
for n„ > 10^^ cm~^ at T = 10^ K jT3]. In general, the higher the temperature and the neutron number 
density are, the better (n, 7) ^ (7, n) equilibrium holds. 
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1.2.2 steady flow and steady /3-flow 

In a steady flow, Y{Z,A) = for all the nuclei in the reaction network and the yields of these nuclei are 
determined by a system of linear algebraic equations. Sample steady-flow r-process calculations can be 
found in p/^. The steady- flow solution also satisfies 

\piZ-l)Y{Z-l) = \p{Z)Y{Z), (8) 

where Y{Z) = Y.aY{Z,A) and \p{Z) = E^Xf^iZ, A)Y{Z, A)/Y{Z) with A^(Z,A) being the total (3- 
decay rate of the nucleus {Z, A). A special case of the steady flow arises when the r-process is also in 
{n, 7) ^ (7, n) equilibrium. In this case, essentially all the abundance in an isotopic chain is carried by 
one or two isotopes with S2n ~ 25*° (see Sec. I1.2.1|l . Such isotopes are called the waiting-point nuclei 
as the r-process must wait for their (3 decay in order to produce the nuclei in the next isotopic chain. 
The above special case of the steady flow is called the steady /5-flow. As can be seen from Eq. (jH}, the 
yield of a waiting-point nucleus in a steady /9-flow is inversely proportional to its /9-decay rate. 

A true steady flow requires that nuclei be fed into the reaction network from below for a sufficiently 
long time. The results from a steady-flow calculation are valid only for certain regions of the network 
where this condition is met. It was found that a steady /3-flow may be obtained for the waiting-point 
nuclei between two magic neutron numbers, e.g., those with 50 < < 82 ^]. If steady /3-flows are 
realized in the r-process, then the peaks at A = 130 and 195 in the solar r-pattern can be attributed to 
the extremely small /3-decay rates of the waiting-point nuclei with the magic neutron numbers A^ = 82 
and 126, respectively. 

1.2.3 fission cycling 

If the r-process involves progenitor nuclei with Z > 80, then fission should be included in the reaction 
network. This may lead to fission cycling, a scenario where the heaviest nucleus produced by the r- 
process fissions and a cyclic flow occurs between this nucleus and its fission fragments in the presence 
of a large neutron abundance. Consider a simple example where (n,7) ^ {1,^) equilibrium holds and 
fission occurs upon the (3 decay of the heaviest waiting-point nucleus with proton number Zf, which 
produces two fragments with fixed proton numbers Zi and Z2 {Z2 > Zi), respectively. A sufficiently 
long time after the onset of the r-process, the nuclei with Z < Zi are depleted to negligible abundances 
by neutron capture and j3 decay. Then the yields of the waiting-point nuclei involved in fission cycling 
are governed by 

YiZ) = \p{Z - l)YiZ - 1) - \p{Z)Y{Z) (9) 
for Zi < Z < Z2 and Z2 < Z < Zf, and 

YiZ,) = \p{Zj)Y{Zf)-\p{Z{)Y{Zr), (10) 
Y{Z2) = MZf)Y{Zf)+\p{Z2-l)Y{Z2-l)-MZ2)Y{Z2). (11) 

With feeding of the nuclear flow at the lower proton numbers Zi and Z2 by fission, a steady state can be 
obtained after a few fission cycles. This results in a rather robust yield pattern with peaks at A ~ 130 
and 195 0. 

1.2.4 freeze-out 

In an r-process event, the neutron number density and the temperature decrease with time. Neutron 
capture ceases to be efficient and the r-process freezes out when 

nn{van,^)r ~ 1, (12) 
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where r is the timescale over which n„ decreases significantly. The freeze-out may involve several stages 
during which the approximations discussed above are most likely to break down. Here full network 
calculations are especially important. After the freeze-out, the progenitor nuclei successively /5-decay 
towards stability. In particular, the progenitor nuclei with magic neutron numbers, which are in the 
peaks of the r-process yield pattern at the freeze-out, /?-decay to the stable nuclei with approximately 
the same mass numbers but with smaller nonmagic neutron numbers. The freeze-out pattern may differ 
significantly from the final yield pattern due to modifications by processes such as /3-delayed neutron 
emission after the freeze-out. 

2 Astrophysical Conditions for the r-Process 

The photo-disintegration rate \^^n{Z,A + 1) is related to the neutron-capture rate n„(f (T„ ,^(2', A)) 
through detailed balance, which leads to Eq. (jH)) for (n, 7) ^ (7, n) equilibrium. The specification of 
the rates for neutron capture and photo-disintegration in the r-process network requires the neutron 
number density n„ and the temperature T. In traditional calculations (e.g., [121 CH CHI CH]), a seed 
nucleus, usually ^^Fe, is irradiated with neutrons at constant values of n„ and T for a time t^j-j.. The 
final yield pattern is then obtained by assuming an instantaneous freeze-out and taking into account the 
modifications of the freeze-out pattern due to /3-delayed neutron emission. Such calculations showed 
that the solar r-pattern must be accounted for by a superposition of the yield patterns resulting from 
distinct sets of n„, T, and tin-. These sets of parameters with typical values of n„ > 10^° cm~^, T ~ 10^ 
K, and > 1 s were considered to represent the astrophysical conditions in the r-process events that 
occurred over the Galactic history prior to the formation of the solar system. 

In a realistic r-process event, n„ and T decrease with time. The new feature of an r-process 
calculation taking this into account can be illustrated by considering the phase during which (n, 7) ^ 
(7, n) equilibrium holds. Due to the evolution of n„ and T, the waiting-point nucleus with mass number 
in the isotopic chain of Z changes with time t, which introduces a time dependence into the effective 
/?-decay rate A^(Z, t) = J2a ^p{Z, A)Y{Z, A)/Y{Z) ^ \fs[Z, Awp(t)]. As A^[Z, Awp(t)] tends to be very 
large at early times when high values of ?t,„ favor the waiting-point nuclei close to the neutron-drip 
line, the duration of an r-process can be greatly reduced. By contrast, the waiting-point nuclei and 
the effective /3-decay rates do not change in the traditional r-process calculations when (n, 7) ^ (7, n) 
equilibrium is assumed. Thus, the neutron irradiation time ti„ used in the traditional calculations to 
obtain the solar r-pattern may be much longer than the actual timescale for an r-process. It is also 
clear that the decrease of Un and T in an r-process event causes (n, 7) ^ (7, n) equilibrium to break 
down eventually and makes the freeze-out more like a multi-stage process rather than a sudden one. 

The evolution of n„ in an r-process event can be described in terms of the neutron-to-seed ratio 

Rn/sit) = r„(t)/n(to), (13) 

where Yn{t) is the neutron abundance at time t and Ys{to) is the initial abundance of seed nuclei. If 
fission can be ignored, conservation of mass and the total number of nuclei gives 

Yn{t)+ J2 ^YiZ,A) = Y„ito) + AsYs{to), (14) 

(Z,A) 

Y{Z,A) = Y,{to), (15) 

{Z,A) 

where As is the mass number of the seed nucleus and the sums extend over all the nuclei with A > Ag. 
Equations (fT^ and (fT3|) can be rewritten as 

A{t) + Rn/sit) =As + Rn/sito), (16) 
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where A{t) = J2{z,a) ^^{Z, A)/ J2(z,a) Y{Z, A) is the average mass number of the nuclei with A> Ag 
that are in the r-process network at time t. In general, when the r-process freezes out aX t = t-poi 
Rn/sitpo) ~ 1 and A(tpo) ^ As + Rn/sito). Thus, the outcome of an r-process event can be simply- 
estimated from the initial neutron-to-seed ratio and the mass number of the seed nucleus. 

The initial neutron-to-seed ratio Rn/s{to) not only provides a convenient means to characterize an 
r-process event but also highlights two important issues: where do the seed nuclei come from and how 
is Rn/sito) determined? Observations of metal-poor stars to be discussed in Sec. IH showed that the 
r-process already occurred in the early history of the Galaxy. This suggests that an r-process event 
cannot rely on some previous astrophysical events to provide the seed nuclei and must produce the 
seed nuclei within the event itself. The production of the seed nuclei and the determination of Rn/s{to) 
in an r-process event can be illustrated by considering a rather generic scenario in which neutron-rich 
material adiabatically expands from high temperature and density. The parameters characterizing the 
expansion that results in a successful r-process can be taken as a new measure of the astrophysical 
conditions for the r-process. 



2.1 Adiabatic expansion from high temperature and density 

The astrophysical input for a nucleosynthesis calculation includes the initial state and the time evolution 
of temperature and density. For electrically-neutral matter expanding from T>10^° K, the initial nuclear 
composition is 

YM = Fe, (17) 
F„(0) = l-n, (18) 

where is the net electron number per baryon, or the electron fraction. Neutron-rich matter has 
Ye < 0.5. For adiabatic expansion, the thermodynamic evolution of the matter is governed by 

TdS + J2 A^(^' A)dY{Z, A) = 0, (19) 

(Z,A) 

where S is the total entropy per baryon for all the particles in the matter (nucleons, nuclei, e~, e+, 
and photons) and fi{Z,A) is the chemical potential (including the rest mass) of the nucleus {Z,A). At 
T > 6 X 10^ K (usually also corresponding to high density), the forward strong and electromagnetic 
reactions are balanced by their reverse reactions. This results in nuclear statistical equilibrium (NSE), 
for which the second term on the left-hand side of Eq. (fT^ vanishes and the adiabatic condition reduces 
to 5" being constant. After NSE breaks down, some entropy is generated by out-of-equilibrium reactions. 
However, this has little effect on considerations of the r-process jTH]. So S may be taken as constant 
during adiabatic expansion. 

If photons are much more abundant than nucleons and nuclei (radiation dominance, 5* ^ 10 in units 
of k per baryon), 

for T>5 X 10^ K when e~ and are also abundant. In Eq. ()20|) . p is the matter density. The coefficient 
3.34 in Eq. (j^Uj) should be replaced by smaller values for T < 5 x 10^ K due to the annihilation of e~ 
and (e.g., the proper coefficient is 1.21 for T < 10^ K). On the other hand, if nucleons and nuclei are 
much more abundant than photons (matter dominance, 5'<10), 5* depends on T and p as ln(T^/^/p). In 
general, S provides a relation between T and p during adiabatic expansion. So only the time evolution 
of T or p needs to be specified explicitly for a nucleosynthesis calculation. One of the parametrizations 
used in the literature ([20], but see [121 HI] for different choices) is 

T{t) = T(0) exp (-t/Tdyn) , (21) 
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where Tdyn is the dynamic timescale of the expansion. Thus, a nucleosynthesis calculation can be done 

dyn 



for matter adiabatically expanding from e.g., T ~ 10^° K once the parameters Y^, S, and r^yn are given. 



2.2 Statistical equilibrium and nucleosynthesis 

In NSE, the abundance of the nucleus {Z, A) can be obtained from fi{Z, A) = Zfip + {A — Z)fin as 



(22) 



where B{Z,A) is the nuclear binding energy of the nucleus {Z,A). The nuclear composition in NSE 
does not depend on the reaction rates and is completely specified by Y^, T, and p together with the 
constraints of charge and mass conservation: 

Ye = Yp+J2 ZY{Z,A), (23) 

iZ,A) 

1 = Y^ + Yp+Yl AY{Z,A), (24) 

where the sums over nuclei exclude the nucleons. After NSE breaks down at T ~ 6 x 10^ K, the nuclear 
composition is governed by the development of quasi-equilibrium (QSE) clusters. The total abundance 
in each QSE cluster is specified by the rates for the nuclear flows into and out of the cluster, but the 
relative abundances of the nuclei within the cluster are the same as given by Eq. ()22|) and again do 
not depend on the reaction rates fI7\ \'23\. A familiar example of QSE is (n, 7) ^ equilibrium, 
for which each isotopic chain is a QSE cluster. The total abundance Y{Z) in the isotopic chain of Z is 
determined by the /3 decay of the nuclei in the isotopic chains of Z — 1 and Z, but the abundance ratio 
Y{Z, A + 1)/Y{Z, A) as given in Eq. © can be simply obtained from Eq. ()22|). 

At 6 X 10^ K > T > 4 X 10^ K, there are two major QSE clusters: the light one involving neutrons, 
protons, and a-particles, and the heavy one involving ^^C and heavier nuclei j2S|- The abundance of 
a-particles in the light cluster is given by Eq. (j22p while the relative abundances of the nuclei within the 
heavy cluster can be estimated by using the same equation. The total abundance in the heavy cluster 
is governed by the net flow out of the light one, which crucially depends on the reactions 3a —>■ + 
a + a + n —>■ ^Be + 7 followed by a + ^Be — ^ ^^C + n, and all the reverse reactions. The nuclear 
composition of the heavy cluster freezes out at T ~ 4 x 10^ K. This composition is slightly modified 
by capture of neutrons and a-particles down to T ~ 2 x 10^ K, at which temperature essentially all 
charged-particle reactions freeze out due to the Coulomb barrier. If there is a sufficiently high neutron 
abundance at this point, (n, 7) ^ (7,'"') equilibrium takes over and an r-process starts. 



2.3 Determination of the initial neutron-to-seed ratio 

As neutron-rich matter with 1^(0) = Ye< 0.5 and Yn{0) = 1— Ye adiabatically expands from T ~ 10^° K, 
its nuclear composition evolves as follows based on the discussion in Sec. 12.21 The first phase of evolution 
is in NSE, during which essentially all the protons are assembled into a-particles. When NSE breaks 
down at T ~ 6 X 10^ K, the nuclear composition is dominated by neutrons and a-particles. During the 
subsequent phase of evolution in QSE, heavy nuclei are produced subject to the infiuence of the bottle- 
neck imposed by the three-body reactions 3a —>■ ^^C + 7 and a + a + n ^Be-l-7 that connect the light 
QSE cluster and the heavy one. The overall neutron excess rj = 1 — 2Ye is crucial to the composition 
of the heavy QSE cluster. For < rj < 0.05 that is similar to the neutron excess r]{Z, A) = 1 — {2Z/A) 
of the individual nuclei in the Fe group, these nuclei are the heavy nuclei favored by QSE as they have 
the largest nuclear binding energy per nucleon B{Z, A)/ A. By contrast, QSE favors the heavy nuclei 
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with Z ~ 35-40 and A ~ 90 (close to the magic neutron number = 50) for ri > 0.05. These nuclei 
have much larger values of ri{Z, A) but similar values of B{Z, A)/ A compared with the Fe group nuclei. 
The overall r] required for producing the heavy nuclei with Z ~ 35-40 and A ~ 90 can be significantly 
smaller than the values of ri{Z, A) ~ 0.1-0.2 for these nuclei if the abundance of a-particles is large. 
This is because rj ~ X„ + Y,{z,A)'Ui.Zi ^) with X„ and X{Z,A) being the mass fractions of 

neutrons and nuclei, respectively. The a-particles have no neutron excess and their presence allows t] 
to be small by reducing the mass fraction of the heavy nuclei that have significant neutron excess. The 
composition of the heavy nuclei freezes out of QSE at T ~ 4 x 10^ K and is modified somewhat by 
capture of neutrons and a-particles down to T ~ 2 X 10^ K. Thereafter, an r-process occurs as the 
heavy nuclei with Z ~ 35-40 and A ~ 90 become the seed nuclei to capture the remaining neutrons. 
The production of heavy nuclei in the presence of significant neutron and a-particle abundances was 
first proposed as the precursor to the r-process by Woosley and Hoffman i24|, who referred to this as 
the "a-process." A thermodynamic approach to understand this process was given in |23j . 

For clarity in discussing adiabatic expansion of neutron-rich matter, the beginning of the expansion 
is denoted as t = 0, that of the a-process is denoted as t^, and that of the r-process is denoted as 
to- The initial neutron-to-seed ratio Rn/sito) for the r-process is determined by the abundances of 
neutrons and heavy nuclei at the end of the a-process, which in turn depend on the net flow into the 
heavy QSE cluster from the light one during the a-process. As discussed in Sec. 12.21 this flow starts 
with the reactions building ^^C. The physics involved in the determination of Rn/s{to) is reflected by 
the dependence of Rn/sito) on Fe, and Xdyn, and can be illustrated by focusing on the reactions 
a + a + n ^ ^Be -|- 7 and a + ^Be ^^C -|- n. For the most part of the a-process, the abundance Yg of 
^Be is approximately given by Eq. ()22j) . As Yp is small during the a-process, it is more convenient to 
use Ya [also given by Eq. (jUJ] to estimate Fg as: 

- ^ ---^^ (w^J (^)'- (^^^) ■ 

The validity of Eq. can be understood from the fragility of ^Be. It only requires a photon of energy 
1.573 MeV for 7 -|- ^Be — a + a + n to occur. Over the temperature range of the a-process, there are 
always some photons on the high-energy tail of the Bose-Einstein distribution that are able to maintain 
the equilibrium between disintegration and formation of ^Be. Similarly, the net flow into the heavy 
QSE cluster is quite small early in the a-process when the temperature is sufficiently high for many 
endothermic reverse reactions such as n + ^^C — > ^Be + a to occur at significant rates. The impediment 
of the fiow to the heavy nuclei is more severe for higher values of S corresponding to more photons per 
baryon. 

The effect of 5* on the a-process becomes even more obvious after the assemblage of neutrons and 
a-particles into heavy nuclei begins in earnest. The rate of assemblage is controlled by the rate for 
the reaction a -|- ^Be — > ^^C + n, which is YgY^^p/ mu) {vaa^ni^Be)) oc [see Eq. (j2Sl)]- For radiation- 
dominated expansion over a given temperature range, this rate is proportional to S^^ [see Eq. (^Uj) ]. 
The integrated assemblage over the duration of the a-process is then proportional to rdyn/5'^. Thus, 
smaller values of Xdyn or higher values of S lead to less consumption of neutrons and less production 
of seed nuclei, both of which tend to increase Rn/sito)- The main dependence of Rn/sita) on Ye is 
through the specification of Yn{ta) at the beginning of the a-process. The nuclear composition at ta is 
dominated by neutrons and a-particles with Ya{toi) ~ Ye/2 and Yn{to) ~ 1 — 2Ye. Clearly, lower values 
of Ye give higher values of Yn{ta)., which means more neutrons at the end of the a-process, and hence, 
higher values of Rn/s(to)- 

The combinations of Yg, S, and Tdyn that give rise to Rn/siio) ~ 100 in matter adiabatically expanding 
from T ~ 10^° K were calculated in [20] and are shown in Fig. 3. For some sets of S and Tdyn, there 
are two possible values of Y^ with one being close to 0.5. This can be explained by the effects of Y^ on 
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Figure 3: Combinations of Y^, S, and Tdyn giving rise to an initial neutron-to-seed ratio of 
Rri/sito) ~ 100 for the r-process in adiabatically expanding matter. Production of nuclei 
with A ~ 195 is expected. See [20] for details. 



the neutron abundance and on the rate for building heavy nuclei during the a-process. A low neutron 
abundance for Y^. ~ 0.5 gives a low rate for building heavy nuclei due to the low abundance of ^Be 
[see Eq. (I25j) ]. This results in a low abundance of the seed nuclei that can compensate for the low 
neutron abundance at the end of the a-process to give Rn/s{to) ~ 100. For a fixed Y^, the required 

1 /3 

values of S and Tdyn to give Rn/s(to) ~ 100 approximately follow the scaling S oc r^y^. This is because 
the integrated assemblage of seed nuclei during the a-process depends on rdyn/5'^. Consequently, for 
the smallest value of Tdyn = 0.0039 s shown in Fig. 3, Rn/s{to) ~ 100 can be obtained over a relatively 
broad range of Ye for S* < 100. The curves for the three very different values of Tdyn shown in Fig. 3 
appear to converge in the region of S* < 10 and Ye < 0.2. This is because for such small values of S 
corresponding to high p, the bottle-neck imposed by the three-body reactions in the production of the 
seed nuclei is no longer important [see Eq. ((23)] • The two major candidate astrophysical environments 
for the r-process to be discussed in Sec. IHl are characterized by 5* ~ 100 (neutrino-driven wind from 
core-collapse supernovae) and S* < 10 (ejecta from neutron star mergers), respectively. 

The combinations of Ye, S, and Tdyn shown in Fig. 3 are considered to represent the conditions for 
an r-process that can produce the nuclei with A ~ 195 from the seed nuclei with Ag ~ 90. Similar 
results were also obtained in [ini ED 12^] • The combinations of Ye and S required to produce the nuclei 
with A ~ 130 would lie to the left of the curves shown in Fig. 3 for fixed values of Tdyn- Illustrative 
r-process calculations with different combinations of Ye, S, and Tdyn were carried out in |Tni I2I| and 

EH]. 
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Figure 4: Neutron number density function of time t for adiabatic expansion with 

Ye = 0.45, Tdyn ~ 0.05 s, and 5* = 130-350. A superposition of these conditions with 
approximately equal weights can account for the solar r-pattern from the peak at A = 130 
up to the peak at A = 195. Note that the dynamic timescale Tdyn used here corresponds to 
a parametrization of T{t) that is different from Eq. (PT|). See [21] for details. 

2.4 Comments on parametric studies 

Both the traditional calculations with constant and T and the more recent calculations based on 
adiabatic expansion of neutron-rich matter from high temperature and density are parametric studies of 
the astrophysical conditions for the r-process. These studies are largely independent of specific models 
and are useful in assessing the potential of any astrophysical environment to be an r-process site. 
An important difference between the traditional calculations and the calculations based on adiabatic 
expansion lies in the production of the nuclei with A < 90-110. These nuclei are produced from the seed 
nucleus ^^Fe by the r-process in the traditional calculations. However, the seed nuclei for the calculations 
based on adiabatic expansion are produced by the a-process prior to the r-process and usually already 
have A ~ 90 (see Sec. 12. 3|) . In fact, significant production by the a-process extends to the nuclei 
with A ~ 110 [211 . Thus, r-process production truly occurs for A > 110 in the calculations based on 
adiabatic expansion. These calculations also shed some new light on the results from the traditional 
calculations. As discussed earlier, the waiting-point nucleus in a given isotopic chain changes due to the 
decrease of n„ and T during the r-process, which means that the parameter t[j.j. used in the traditional 
calculations does not represent the actual timescale for an r-process. On the other hand, adiabatic 
expansion that can produce the peaks at A = 130 and 195 in the solar r-pattern was shown to result in 
steady /3-fiow at the freeze-out of the r-process with values of n„ and T similar to those obtained from 
the traditional calculations . Sample time evolution of n„ for adiabatic expansion with fixed Yg and 
Tdyn but different values of S [21] is shown in Fig. 4. 

Some comments on the use of the solar r-pattern in parametric studies are in order. The traditional 
calculations require that a suitable superposition of the yield patterns resulting from different combi- 
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nations of ?t,„, T, and ti„ match the solar r-pattern shown in Fig. 2 plus a peak at A = 80. There 
are several potential problems with this requirement. First, the solar "r-process" abundances of the 
nuclei above the Fe group but with A < 90-110 were derived by subtracting the calculated s-process 
contributions from the total solar abundances jH] and may contain substantial contributions from the 
a-process, which can produce many of these nuclei [21]. Therefore, the results based on fitting the solar 
"r-pattern" associated with the peak at A = 80 and possibly up to the nuclei with A ~ 110 must be 
viewed with caution. Furthermore, as will be discussed in Sec. EJ the basic templates underlying the 
overall solar r-pattern may be more complicated than those employed in the traditional calculations. 
The r-process yield pattern also depends on the nuclear physical input, almost all of which are calcu- 
lated from theory with little experimental guidance. In view of the above issues, parametric studies may 
be of more general use if they are based on the simple criterion that an r-process must have an initial 
neutron-to-seed ratio Rn/siio) = {A) — Ag in order to produce nuclei with an average mass number {A) 
from the seed nuclei with mass number Ag. In contrast to the case of the r-process, the essential nuclear 
physical input is known for the a-process that determines Rn/s{to) and Ag. 

Parametric studies of the r-process have the advantage that no specific astrophysical models are 
required. This also leads to some serious flaws. For example, the parametric studies based on adiabatic 
expansion cannot give the absolute amount of material that has a specific combination of Y^,, S, and 
Tdyn- Furthermore, a real r-process environment may be improperly characterized by some parametric 
studies. This flaw has already been demonstrated by the traditional calculations with constant n„ and 
T. In addition to the problem with the parameter ti„, the freeze-out of the r-process is treated as 
instantaneous in such calculations. However, the actual freeze-out is a multi-stage process during which 
the r-process yield pattern may be smoothed independent of /3-delayed neutron emission (e.g., pTll27j ). 
It was also proposed ^29^ that effects during the freeze-out could be responsible for the small peak in 
the rare earth region {A ~ 160) of the solar r-pattern (see Fig. 2). Such details of the freeze-out are 
revealed only when the time evolution of n„ and T in the r-process event is taken into account. In many 
respects, the calculations based on adiabatic expansion are more realistic. However, it is quite possible 
that they still miss some important ingredient. For example, additional parameters characterizing the 
neutrino fluxes may be required if the r-process occurs in core-collapse supernovae (see Sec. I3.1.H) . 

3 Astrophysical Models of the r-Process 

A neutron-rich environment is required for an r-process. It is interesting to survey the possibilities to 
obtain such environments during the evolution of the universe. After the big bang, the universe goes 
through some phase transitions, which might produce density inhomogeneities of such large sizes that 
only neutrons but not protons can diffuse to uniformity prior to the onset of big bang nucleosynthesis. 
This segregation of neutrons and protons could result in neutron-rich regions where a primordial r- 
process might occur (e.g., jHOJ EII)- The details of such a scenario will not be discussed here as its 
prediction of a minimum r-process enrichment lacks observational support. The major outcome of big 
bang nucleosynthesis is an overall proton-rich nuclear composition with ^ 76% of H and ~ 24% of ^He 
by mass. This composition has hardly changed over the history of the universe. Subsequent to big 
bang nucleosynthesis, neutron excess can be produced locally by processes of weak interaction inside 
stars. Indeed, the small neutron excess generated by (3~^ decay during evolution of stars with masses of 
M > 10 Mq is crucial to the production of the solar abundance pattern for the nuclei from ^^Na to the 
Fe group (e.g., jS2])- Weak interaction also participates in the reaction sequence of p + ^^C ^^N + 7, 
^^N ^^C + e"*" + i^e, and a + ^^C —>■ + which is the main neutron source for the s-process in AGB 
stars (e.g., IH]). The most dramatic production of neutron excess occurs during the deleptonization of a 
protoneutron star that is produced by the gravitational collapse of a stellar core. Over a period of ~ 10 
s, the protoneutron star loses nearly all of its initial electrons and protons through the neutronization 
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reaction e~ + p ^ n + i'e- Concurrently with the deleptonization, the protoneutron star also releases its 
enormous gravitational binding energy through emission of i/g, Ug, u^, z/^, u^-, and Ur- The interaction 
between the neutrinos and the material above the protoneutron star drives a supernova explosion on a 
timescale of < 1 s (e.g., |33j) and ejects a small amount of material in a wind that lasts for a period 
of ~ 10 s (e.g., [33]). The neutron- richness of the wind material is set by the competition between the 
reactions Ug + n — > p+e~ and I'e+P n + e~^ (e-g-, P31 l36j). The high neutron number density of > 10^° 
cm~^ required for the r-process (see Sec. |2] and Fig. 4) may be obtained in the neutrino-driven wind 
associated with protoneutron star formation in core-collapse supernovae. Material with such neutron 
number density may also be ejected when an old neutron star is disrupted during its merging with 
another neutron star or a black hole. The core-collapse supernova and neutron star merger models of 
the r-process are discussed in detail below. 

3.1 Core-collapse supernova models of the r-process 

The evolution, explosion, and nucleosynthesis of massive stars are a rich subject (see jSTj for a recent 
review) . Only a simple sketch of the evolution and explosion of such stars is given here to provide some 
general astrophysical context for the discussion of the r-process. Through various stages of hydrostatic 
burning, massive stars with M > 10 Mq develop an onion-skin structure with an Fe core surrounded by 
shells of successively lighter elements from Si to H. As no more nuclear binding energy can be released by 
burning Fe, the Fe core undergoes gravitational collapse. When the inner core reaches nuclear density, 
it becomes a protoneutron star and bounces. This launches a shock, which fails to exit the outer core 
mainly due to energy loss from dissociation of Fe. Although a successful mechanism remains to be 
demonstrated (see [3Hj for a review), it is considered that the shock is revived by the neutrinos emitted 
from the protoneutron star and proceeds to make a Type II supernova (SN II). Massive stars with 
8 Mq < M < 10 Mq develop a degenerate 0-Ne-Mg core with insignificant shells jnS] • The core collapse 
in this case is triggered by electron capture jlUj and again leads to protoneutron star formation and an 
SN II |¥T1 142j . A bare white dwarf in a binary may also collapse into a protoneutron star due to mass 
accretion from its binary companion |43ij. The supernovae resulting from accretion-induced collapse of 
white dwarfs are sometimes referred to as the "silent supernovae" due to the lack of prominent optical 
display. 

Whether formed by Fe core collapse, 0-Ne-Mg core collapse, or accretion-induced collapse (AIC) of a 
white dwarf, the protoneutron star releases its gravitational binding energy through neutrino emission, 
the characteristics of which have been modeled by numerical transport calculations (e.g., |^-[49]). 
The gravitational binding energy of a protoneutron star with a mass of Mns ~ 1.4 and a radius of 
i?NS ~ 10 km is ~ (3/5)(j'M^g/i?NS ~ 3 x 10^^ erg, where G is the gravitational constant. This energy 
is emitted approximately equally in each neutrino species over the neutrino diffusion timescale. All 
neutrino species have neutral-current scattering on nucleons and diffuse out of the protoneutron star 
on a timescale of ~ 10 s. The neutrino luminosities satisfy 

L^, ~ Lp, ~ L^^ ^ Lp^ ^ L^^ ^ Lp^ (26) 

and have a typical value of ~ 10^^ erg s~^ for each species. During diffusion, different neutrino species 
stop exchanging energy with matter in different decoupling regions, the temperatures of which charac- 
terize the neutrino energy spectra emergent from the protoneutron star. Energy exchange with matter 
occurs mainly through neutral-current scattering on electrons for u^, z/^, u^-, and Ur, but mainly through 
the charged-current reactions u^. + n ^ p + e~ and Uf. + p —>■ n + e'^ with higher efficiency for i/g and z/g. 
In fact, energy exchange is more efficient through Ue + n ^ p + e~ than through + p n + as 
there are more neutrons than protons inside the protoneutron star. Consequently, the decoupling from 
local thermodynamic equilibrium with matter occurs first for z/^, z/^, Ur, and Ur, next for z/g, and last 
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for i/g. Thus, the average energies corresponding to the emergent neutrino spectra satisfy 



{EuJ ^ {E,^) ^ {E..) ^ {E,^) > {E,J > {E,J (27) 

with typical values of {E^J ^ 11 MeV, {EpJ ^ 16 MeV, and (E^^) ^ 25 MeV. 

Close to the protoneutron star, the temperature is T > 10^° K and the material is in NSE with its 
nuclear composition dominated by nucleons. This material is heated by the reactions z/g + n — > p + e~ 
and + p —>■ n + e'^ and expands away from the protoneutron star as a neutrino-driven wind. Above 
a radius corresponding to T ~ 6 x 10^ K, the neutrino reactions become inefficient due to the ever- 
decreasing neutrino flux. The wind material then essentially undergoes adiabatic expansion with flxed 
values of Ye, S, and Tdyn as in the parametric r-process study described in Sec. |21 The parameters Ye, 
S, and Tdyn in the wind are determined by the history of neutrino interaction with the wind material 
at T > 10^° K [36 . In particular. Ye is set at sufficiently low values of T so that the rates Xe+n and Ae-p 
for the reactions + n p + Ue and e~ + p ^ n + Ue are unimportant compared with the rates X^^n 
and Xp^p for the reactions Ue + n p + e~ and Ue + p —>■ n + . This gives an estimate of Ye as 



Yg ^ ^ = - . (28) 

Xu^n + Xp^p 1 + (Ap^p/ X^^n) 

With L^^ ~ Lp^ and {EpJ > {E^J, it was shown that Ap^p > X^^n, and hence. Ye < 0.5 can be obtained 
at least for some part of the period of ~ 10 s over which the neutrino-driven wind occurs jHE]- This 
wind was first proposed as a site of the r-process in |50j. However, only one SN II model gave the 
adequate values of Ye, S, and Tdyn for producing the peaks a.t A = 130 and 195 in the solar r-pattern 
jni]. The typical values of 0.45 <Ye< 0.5, S ~ 100, and 0.01 s < Tdyn ~ 0.1 s obtained in another SN 
II model j2ni and in analytic and numerical studies of the neutrino-driven wind j36l 1221 were unable to 
give an r-process (see Fig. 3). Parametric studies of r-process nucleosynthesis in the neutrino-driven 
wind were carried out in 



3.1.1 effects of neutrinos 

As will be discussed in Sec. 14.31 in order to account for the solar r-process abundances associated with 
the peaks at A = 130 and 195, each supernova must eject ~ 10^^-10"^ Mq of r-process material. 
Although the current neutrino-driven wind models have difficulty in providing the r-process conditions, 
the wind naturally ejects ~ 10~^-10~^ Mq of material over a period of ~ 1 s. This is because the 
small heating rate due to the weakness of neutrino interaction permits material to escape from the deep 
gravitational potential of the protoneutron star at a typical rate of ~ 10~^-10~^ Mq jnnilS2]- Indeed, 
the ability to eject a tiny but interesting amount of material was recognized as an attractive feature of 
the neutrino-driven wind model of the r-process (e.g., [26j). On average, a nucleon obtains ~ 20 MeV 
from each interaction with Ug or z/g. In order to escape from the protoneutron star gravitational potential 
of GMns'^wZ-Rns ~ 200 MeV, a nucleon in the wind must interact with z/g and Ue for ~ 10 times. This 
suggests that the effects of neutrino interaction may be important even during the essentially adiabatic 
expansion of the wind material. Some of these effects are discussed below in connection with r-process 
nucleosynthesis. 

At T > lO^'' K, the evolution of Ye in the wind material is governed by 

Ye = {X^^n + Ae+n)^n ~ {Xp^p + Ag-p)!^ = X^^n + Ae+„ — {X^^n + Ap^p + Ag+„ -|- Ae-p)yg, (29) 

where y„ = 1 — Fg and Yp = Ye have been used to obtain the second equality. For small values of Tdyn, 
Yg freezes out with a value given by Eq. (j2Hl) when T is so low that Ag+„ and Ag-p can be neglected 
compared with Xu^n and Ap^p but nucleons still dominate the nuclear composition. However, if neutrino 
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interaction continues to be significant when NSE starts to favor a-particles over nucleons at T < 10 
K, tlie equation governing tlie evolution of Yf, becomes 




(30) 



which is obtained by neglecting Ae+n and Ae-p and using Yn = I — Ye — (Xq,/2) and Yp = Ye — {Xa/2), 
with Xa being the mass fraction of a-particles, in the first equahty of Eq. (j29|) . In the presence of 
a-particles, the evolution of Ye tries to reach 



which is increased from the value in Eq. for Ap^p > X^^n- This so-called a-effect was first discussed 
in |53]. Depending on Tdyn, the a-effect may significantly decrease the initial neutron-to-seed ratio and 
hinder the r-process [21] • The effects of neutrino interaction with nucleons and nuclei on the evolution 
of Ye were studied in detail in (HH! ■ 

For large values of S and Xdyn, the a-particles present an additional problem due to their own 
interaction with neutrinos. In particular, z/^, i/^, i/^, and z/^ with the highest average energy can induce 
proton spallation on the a-particles. The daughter nucleus from this process provides a new path 
to produce the seed nuclei starting with the reactions a -|- ^ ''Li -|- 7 and a + "^Li ^^B + 7. As 
this path does not involve the inefficient three-body reactions to burn the a-particles (see Sec. 12. 2|) . the 
production of the seed nuclei can be significantly enhanced, which reduces the initial neutron-to-seed 
ratio and again hinders the r-process 

Provided that the above two problems can be avoided and an r-process occurs in the neutrino-driven 
wind, neutrino interaction may have some direct effects on r-process nucleosynthesis. First of all, the 
progenitor nuclei encountered during the r-process are extremely neutron-rich and have large cross 
sections for z/e capture (see [HZl EHI for recent calculations). In contrast to P decay, Ue capture on the 
progenitor nuclei are insensitive to the presence of magic neutron numbers. Because the abundance 
peaks at A = 130 and 195 in the solar r-pattern are usually attributed to the slow /?-decay rates of the 
progenitor nuclei with the magic neutron numbers = 82 and 126, respectively, no such peaks would be 
produced if the rates for Ue capture were to dominate those for (3 decay at the freeze-out of the r-process 
j53j . However, as the r-process proceeds in the wind material, this material is also expanding away from 
the protoneutron star. Consequently, the neutrino flux, and hence, the neutrino interaction experienced 
by the wind material decreases during the r-process. It is possible that z/e capture is significant in the 
early phase of the r-process but becomes negligible compared with j3 decay at the freeze-out |^ ItiOj . 
In this case, Ue capture can accelerate the progress from one isotopic chain to the next and reduce the 
duration of the r-process as proposed in [61j. The effects of Ue capture on the r-process flow were also 
discussed in [^21 EBI • 

Neutrino interaction may also be important during decay towards stability after the r-process freezes 
out. Both Ue capture and neutral-current reactions with u^, u^, u^-, and z/^ can greatly excite the 
progenitor nuclei, which then deexcite through neutron emission (see EHl ioi recent calculations) 
or possibly fission. It was shown that the solar r-process abundances of the nuclei with A = 124-126 
and A = 183-187 may be completely accounted for by neutrino-induced neutron emission from the 
progenitor nuclei in the abundance peaks at A = 130 and 195, respectively [Sni El]- This neutrino 
postprocessing effect is shown for the mass region near A = 195 in Fig. 5. This effect also constrains 
the total exposure to neutrinos after the freeze-out as too much neutrino postprocessing overproduces 
the nuclei below the abundance peaks at A = 130 and 195 [W, "M]. The effect of neutrino-induced 
fission during decay towards stability will be discussed in Sec. 14.3.21 in connection with the r-patterns 
observed in metal-poor stars. 




(31) 
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Figure 5: Effect of neutrino postprocessing for the mass region near the abundance peak at 
A = 195. The abundances before and after the postprocessing are given by the sohd and 
dashed curves, respectively. The filled circles (some with error bars) give the solar r-process 
abundances. The region with A = 183-187 is highlighted in the inset. See [HHlini] for details. 



With a number of possible effects of neutrinos, it is important to do a self-consistent study that 
includes all the neutrino interaction processes. A significant step in this direction was taken in |54j . 
where a comparison was made between the results calculated by including neutrino interaction to various 
extent during adiabatic expansion of the wind material. However, this comparison was based on a fixed 
set of wind parameters and the neutrino luminosities used to calculate the neutrino reaction rates were 
not compatible with the parameter Xdyn adopted for the wind. In a fully self-consistent study, the wind 
parameters such as 5* and Tdyn should be determined by the same neutrino emission characteristics that 
are used to study the neutrino effects on r-process nucleosynthesis. It is essential that such a study be 
carried out in the future. 

3.1.2 remedies for the neutrino-driven wind model of the r-process 

As will be discussed in Sec. EJ a number of observations are in support of core-collapse supernovae being 
the major site of the r-process. The amount of r-process material required from each supernova to ac- 
count for the solar r-process abundances associated with the peaks at A = 130 and 195 is approximately 
the same as the amount of ejecta in the neutrino-driven wind. However, current models fail to provide 
the conditions for an r-process to occur in the wind. This failure may simply reflect the uncertainties 
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in the models and can be remedied when better physical input is used. A major uncertainty in the 
wind models concerns the properties of hot and dense matter, which affect the protoneutron star mass 
Mns and radius -Rns- General relativistic effects are crucial when 2GMns/(-Rnsc^) (c being the speed of 
light) approaches unity. For the nominal values of Mns = 1-4 Mq and i?NS = 10 km, such effects already 
lead to significantly larger values of S and smaller values of Tdyn than given by Newtonian calculations 
for the neutrino-driven wind Both increase of S and decrease of Tdyn tend to increase the initial 
neutron-to-seed ratio for the r-process (see Sec. I2.3|l . Thus, massive and/or compact protoneutron stars 
with large values of 2GMns/(-Rnsc^) may have a wind with adequate conditions for an r-process (see 
ISSIIH2] and [HSj^inZI)- The properties of hot and dense matter also affect neutrino interaction processes 
in protoneutron stars, and hence, the neutrino emission characteristics obtained from transport calcu- 
lations (e.g., IHHl^IZni)- The parameter Y^. for the neutrino-driven wind is sensitive to the difference 
between the rates X^^n and Ap^p (see Eq. |2H])5 which are determined by the emission characteristics of 
z/g and z/g |35| ESI 1^ • It remains to be seen if low values of Y^. favorable for an r-process might result 
from neutrino transport calculations based on better understanding of the properties of hot and dense 
matter. 

Another major uncertainty in the neutrino-driven wind models concerns the properties of neutrinos. 
As mentioned above. Ye for the wind is sensitive to the difference between the rates X^^n and Ap^p. 
Because the average energy of z/^, z/^, z/,-, and Vr is higher than that of either z/g or z/g [see Eq. (l27j) ]. 
neutrino flavor transformation between the former group and the latter can greatly affect X^^n and Ap^p, 
and hence, Yg jHE]- There might also exist a sterile neutrino that has no standard weak interaction. 
Flavor transformation involving the sterile neutrino can dramatically decrease Yg in the neutrino-driven 
wind by reducing the z/g flux ^72i||73j. With the z/g flux greatly reduced in the region where a-particles 
are being formed, the a-effect is also rendered inoperative [73 US] . Therefore, sterile neutrinos may 
play an essential role in enabling r-process nucleosynthesis in the neutrino-driven wind. Evidence for 
the existence of such neutrinos may come from future experiments on neutrino oscillations. 

In typical neutrino-driven wind models, various neutrino interaction processes provide the energy for 
driving the wind. The radial distribution of these energy sources has significant effects on the parameters 
S and Tdyn for the wind |36j- A broad distribution above the radius at which the total energy deposition 
rate peaks tends to favor an r-process by increasing S and decreasing Tdyn [36, 52 . Such a distribution 
may be obtained when neutrino flavor transformation or especially some extra energy source other than 
neutrino interaction (e.g., magnetic field and rotation of the neutron star [SF) is included in the wind 
models. 

In considering the neutrino-driven wind model of the r-process, it is also important to note the 
difference between the winds associated with Fe core collapse, 0-Ne-Mg core collapse, and AIC of a 
white dwarf. Although the mass ejection rate and S in the wind are determined near the protoneutron 
star, the structure of the wind further away may be affected by the outer boundary imposed by the 
supernova (e.g., [Sn])- For example, the reverse shock in the Fe-core-collapse events may slow down 
the wind significantly pUl whereas no such shock is present in the O-Ne-Mg-core-coUapse or AIC 
events. This may have important consequences for r-process nucleosynthesis associated with core- 
collapse supernovae. 

3.1.3 other core-collapse supernova models of the r-process 

The neutrino-driven wind model of the r-process is closely related to the neutrino-driven mechanism of 
core-collapse supernovae. It is recognized that this mechanism crucially depends on neutrino transport 
and convection (see jSHI for a review), a simultaneous and full treatment of which is very difficult to 
implement numerically. Thus, the robustness of this mechanism remains to be demonstrated. Two 
alternative explosion mechanisms were proposed for core-collapse supernovae and both have associated 
r-process scenarios. In the prompt mechanism, the bounce of the inner core upon reaching nuclear 
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density launches an energetic shock, which then drives a successful explosion (e.g., |SH1)- This mechanism 
might possibly work for the accretion-induced collapse of white dwarfs and the collapse of 0-Ne-Mg 
cores and maybe some low-mass Fe cores. The inner ejecta from a prompt explosion is neutron-rich and 
has been studied as a possible site of the r-process jlTl[711[73]. However, the conditions in this ejecta 
and the amount of ejected r-process material are quite uncertain as the prompt mechanism underlying 
this r-process model is even more problematic than the neutrino-driven mechanism. Another supernova 
mechanism relies on the formation of magnetohydrodynamic jets following the core collapse (e.g., jTHj- 
[ZHj)- While r-process nucleosynthesis associated with jets was discussed quite some time ago [10], an 
intriguing new scenario was proposed by Cameron jHO] recently. In this scenario, jets are associated 
with the formation of an accretion disk around the protoneutron star. The accretion disk has such 
high densities that electrons are degenerate but neutrinos can escape. The chemical equilibrium among 
neutrons, protons, nuclei, and electrons results in large abundance ratios of neutrons to heavy nuclei in 
the disk. The r-process was considered to occur as material is transported from the disk to the base of 
the jets and ejected. Clearly, more quantitative studies of this r-process model are worth pursuing. 

3.2 Neutron star merger models of the r-process 

Two binary neutron star (NS-NS) systems, PSR 1913-1-16 |81j and PSR 1534-1-12 [H21, were observed 
in the Galaxy. The neutron stars in an NS-NS binary eventually merge due to orbital decay caused 
by gravitational radiation. The total time from birth to merger is ~ 4 x 10*^ yr for PSR 1913-1-16 and 
^ 3 X 10^ yr for PSR 1534-1-12 [HH]. Estimates for the rate of NS-NS mergers in the Galaxy range from 
~ 10~^ to ~ 3 X 10~^ yr~^ with the best guess being ~ 10"^ yr~^ (e.g., jHSI^IHni)- The birth rates of 
neutron star-black hole (NS-BH) and NS-NS binaries are comparable. However, the fraction of NS-BH 
binaries having the appropriate orbital periods for merging within the age of the universe (~ 10^° yr) is 
uncertain due to their complicated evolution involving mass exchange [SSj- In any case, the total rate 
of neutron star (including NS-NS and NS-BH) mergers in the Galaxy is perhaps ~ 10~^ yr~^, which 
is ~ 10^ times smaller than the Galactic rate of SNe H [Sn]. This means that each merger must eject 
> 10~^ Mq of r-process material if neutron star mergers were solely responsible for the solar r-process 
abundances associated with the peaks aX A = 130 and 195 (~ 10~^-10~^ Mq of r-process material is 
required from each event in the case of core-collapse supernovae). 

It was estimated that < 5% of the original neutron star mass may be ejected during tidal disruption 
of the neutron star in an NS-BH merger |H71|HH]. Estimates for the amount of cold neutron star matter 
ejected during an NS-NS merger range from ~ 1O~^M0 to ~ 4 x 10~^ [Hn]-|n2]- The ejected cold 
neutron star matter is decompressed and heated to high temperature (>10^ K) by P decay and possibly 
fission jnSl Ell- Further evolution of the decompressed neutron star matter is similar to the adiabatic 
expansion of neutron-rich matter described in Sec. |21 and could lead to an r-process. The r-patterns 
resulting from decompression of cold neutron star matter with = 0.05, 0.10, 0.12, and 0.15 jHll are 
shown in Fig. 6. For these values of Yf,, fission cycling occurs and only the heavy r-process nuclei 
with A > 130 are produced. Neutron star mergers over the Galactic history would make significant 
contributions or even completely account for the solar r-process abundances at A > 130 (e.g., |S71 ISH] 
and [Oni^inS]) if ~ 10"'^ Mq of cold neutron star matter with Fe ~ 0.1 could be ejected from each event. 
The implications of meteoritic data and stellar observations for this r-process model will be discussed 
in Sec. HI 

Several issues regarding the neutron star merger model of the r-process need to be addressed. First 
of all, numerical calculations gave a wide range for the amount of cold neutron star matter ejected 
during an NS-NS merger. It was found that this amount sensitively depends on the initial spins of 
the neutron stars and on the equation of state at high density jHOl^inS- General relativity may also 
be important. All these issues should be examined by future numerical studies with higher resolutions 
(a high- resolution Newtonian calculation was recently carried out in ^^). Simulations of an NS-BH 
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Figure 6: The r-patterns (solid curves) resulting from decompression of cold neutron star 
matter with Yg = 0.05, 0.10, 0.12, and 0.15, which may be ejected from neutron star mergers. 
The solar r-pattern (pluses) is also shown for comparison. See |M| for details. 



merger are also needed to determine the actual amount of ejecta in such events. The other major 
uncertainty in the neutron star merger model of the r-process concerns in the ejecta. The values 
of not only depend on which part of the neutron star is ejected but may also be affected by weak 
interaction with e~, e"^, z/g, and z/g. Although one study indicated that weak interaction has little effect 
on Ye and the ejecta essentially retains the Yg values of the cold neutron star matter, such values are 
still sensitive to the equation of state at high density |22j. In general, as in the case of core-collapse 
supernovae, much work with improved numerical methods and physical input is required to resolve the 
outstanding issues regarding neutron star mergers as a possible site of the r-process. 



3.3 Similarities between core-collapse supernova and neutron star merger 
models of the r-process 

Despite many differences between core-collapse supernovae and neutron star mergers, the neutrino- 
driven wind in the former events and the cold neutron star matter ejected in the latter events follow 
similar thermodynamic evolution, which may lead to the same generic picture for the r-process as 
described in Sec. 12.31 Furthermore, a neutrino-driven wind can also occur in neutron star mergers 
|Hni- While i/g, z/g, z/^, z/^, z/,-, and Ur emitted during such events have approximately the same energy 
spectra as in the case of core-collapse supernovae, the luminosity of z/g is ~ 3 times larger than that 
of z/g due to the dominant production of z/g via the reaction + n ^ p + Ue lEHl 102] • Consequently, 
the neutrino-driven wind from neutron star mergers can have Y^ values much less than 0.5 [see Eq. 
f l28l) ]. It was estimated that ~ lO'^-lO'^ Mq of very neutron-rich material may be ejected in this wind 
||89|j . However, detailed studies of this wind require much better treatment of neutrino transport than 
implemented in the existing models 1^2]. It is important that such studies be carried out in the future. 
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4 Observational Studies of the r-Process 



Enormous progress has been made in observational studies of the r-process since the last major review 
by Cowan et al. jH] in 1991. This includes meteoritic evidence for the diversity of r-process sources and 
demonstration of some regularity in r-patterns and of large dispersions in r-process abundances at low 
metallicities by observations of metal-poor stars. These observational advances and their implications 
for the r-process are discussed in connection with a three-component model for abundances in metal- 
poor stars. Other models addressing the abundances in metal-poor stars and their implications for the 
r-process are also described. 



4.1 Meteoritic data and diverse sources for the r-process 

Meteorites were formed during the birth of the solar system ^ 4.6 x 10^ yr ago. They provide important 
information on the inventory of radioactive nuclei in the early solar system (ESS). For example, ^^^Hf has 
a lifetime of ri82 = 1-3 x 10'' yr and /^-decays to ^^^Ta, which in turn /^-decays with a lifetime of 165 days 
to the stable nucleus ^^^W. Any ^^^Hf that had been incorporated into the meteorites decayed to ^^^W 
long ago. Unlike ^^^W, the stable nuclei ^^°Hf and ^^^W do not receive any contributions from radioactive 
decay. As ^'^^Hf and ^^°Hf are chemically identical, a meteorite with a high abundance of ^^°Hf also 
has a high initial abundance of ^^^Hf. Therefore, if a significant inventory of ^^^Hf existed in the ESS, 
there would be a linear correlation between the present abundance ratios ^^^W/^^^W and ^^''Hf/^^^W 
in the meteorites. Such a correlation was found and gave (^'^^Hf/^^°Hf)Ess ~ 10~^ for the abundance 
ratio of ^^^Hf to ^^°Hf in the ESS [97i ^Hj (see also PU]- |1U2] ). Similarly, meteoritic measurements gave 
(^i29j/i27j)^gg ^ ^Yie abundance ratio of ^^^I (with a lifetime of r^g = 2.3 x 10^ yr) to ^^"^l 

(stable) in the ESS [103, 104J. 

The nucleus ^^^Hf /3-decays with a lifetime of 61 days in the laboratory and cannot provide a 
significant branching for neutron capture to produce ^^^Hf during the s-process. So major production 
of ^^^Hf can only occur in the r-process. Both ^^^I and ^^^I are essentially pure r-process nuclei. The 
meteoritic data on ^^^Hf and ^^^I could not be explained if there were only a unique r-process source 
with a universal yield pattern. This can be shown by the method of contradiction. If ^^^Hf and ^^^I 
were produced with a yield ratio of 1^129/ ^i82 by a unique r-process source, then 
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This is because ^^^I survives longer than ^'^^Hf. The yield ratio 1^129/^182 can be rewritten as 



YlS2 1^127 \ ^^^W. 



(33) 







where the yield ratio ^127/^132 has been set to be the same as the solar r-process abundance ratio of 
^^^I to i^^W (with ^^"^Wr being the part of ^^^W contributed by the r-process) under the assumption of 
a universal r-process yield pattern. Equations (fH^ and can be rearranged into 

r'ns. , (34) 



(182Hf/180Hf)Ess F127 V'^'W, 

The left-hand side of Eq. (I34p is ~ 1 based on the meteoritic data, while the right-hand side is ~ 4 
by taking ^129/^127 = 1-4 and (^^^m/^^^Wr)^ = 0.0541/0.019 HHn]. Thus, Eq. §^ is violated by a 
factor of ~ 4 and the underlying assumption that ^'^^Hf and ^^^I were produced by a unique r-process 
source must be invalid. 
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While the diversity of r-process sources can be estabhshed by using only the meteoritic data on ^^^Hf 
and ^^^I, more may be learned about the nature of these sources by making the following assumption. 
Suppose that ^^^Hf and ^^^I are produced by two distinct kinds of r-process events, which occurred at 
regular intervals of Anf and Aj, respectively, over a uniform production period of Tup ~ 10^° yr prior to 
the formation of the solar system. If the interval between the last event of each kind and the formation 
of the solar system is exactly Anf or Aj, then 




ESS 



ESS 



Tl82 \ Anf/ ri82 



Tup/ exp(AHf/ri82) - 1' 

^129 \ /T129\ Ai/ri29 



Y127J VTup/ exp(Ai/ri29) - 1' 



(35) 
(36) 



Equations and (jHUj) are valid for Anf, Aj <^ Tup, in which case the radioactive nuclei in the 
ESS were contributed by the last few events prior to the formation of the solar system while the 
stable nuclei were accumulated over essentially the entire period of Tup. With (^^^Hf/^^^Wr.)Ess = 
(^S2Hf/i80Hf)Ess(^^°Hf/i«2^^)^ = 2.8 X 10-^ and Yug/Yuj = 1.4, Eqs. ^ and ^ can be solved 
to give Anf = 3.3 x 10^ yr and Ai = 1.2 x 10® yr, which should indeed represent two distinct kinds 
of r-process events. Based on similar arguments, Wasserburg et al. |lU5j first pointed out that the 
meteoritic data on ^^^Hf and ^^^I require diverse sources for the r-process (see also |lU6p . In general, 
other nuclei with A > 130 are produced together with ^^^Hf while those with A < 130 are produced 
together with ^^^I. The meteoritic data on ^^^Hf and ^^^I suggest that there must be at least two distinct 
kinds of r-process events: the high-frequency H events producing mainly the heavy r-process nuclei 
with A > 130 {H stands for "high-frequency" and "heavy r-process nuclei") and the low-frequency 
L events producing dominantly the light r-process nuclei with A < 130 {L stands for "low-frequency" 
and "light r-process nuclei") |6( H I105| I106j . An average interstellar medium (ISM) is enriched with the 
appropriate r-process nuclei at a frequency of fn = A^f ~ (3 x 10'' yr)^^ by the H events and at a 
frequency of /l = Af ^ ~ (10® J^)~^ by the L events. 

The frequencies of the H and L events can be compared with the frequencies for replenishment of 
newly-synthesized material in the ISM by core-collapse supernovae and neutron star mergers. The total 
kinetic energy of the ejecta from a supernova or neutron star merger is typically ~ 10^^ erg. The amount 
of ISM required to dissipate the energy and the momentum of the ejecta is Mmix ~ 3 x 10*^ M0 (e.g., 
|107j ). For a present event rate of /g in the Galaxy corresponding to a total gas mass of Mgas ~ 10^° Mq, 
the frequency for replenishment of newly-synthesized material in the ISM by the relevant events is 



/mix = (10' yr)- 



G 



(30 yr)- 



3 X lO^M^ / V M 




(37) 



mix 



The present Galactic rate of ~ (30 yr)"^ for core-collapse supernovae [SHI corresponds to /^^ 
(10^ jt)~^ while that of ~ (10^ yr)^^ for neutron star mergers (e.g., |HS1) corresponds to f^f^ ~ 
(3 X 10^° yr)^^- As fH>fL^ /mix^) it seems that neutron star mergers cannot be associated with 
either the H or L events. This conclusion is unaffected if the upper limit of ~ (10^ yr)"^ jHll instead 
of the nominal value of ~ (10^ yr)^^ is used for the present Galactic rate of neutron star mergers. By 
contrast, is comparable to fn and fi within a factor of ~ 10, which suggests that the H and 
L events may be associated with core-collapse supernovae |1051 110^ . It is important that the above 
argument be substantiated by more detailed studies on the mixing of the ejecta from core-collapse 
supernovae and neutron star mergers with the ISM and on the occurrences of these events over the 
Galactic history. Independent of the association with specific astrophysical sites, the characteristics 
of the H and L events are further tested and elucidated by observations of abundances in metal-poor 
stars. 
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4.2 Observations of abundances in metal-poor stars 



The effects of diverse sources for the r-process would be most prominent in the early Galaxy where not 
many r-process events had contributed to the abundances in a reference ISM. The chemical evolution 
of the early Galaxy is reflected by the abundances in metal-poor stars that reside in the Galactic halo. 
Observations of abundances in a large number of metal-poor stars as well as detailed studies covering 
many elements in individual stars have been carried out by a number of groups (e.g., |lU9j - [TT^ ). 
These observations provide strong support for the diversity of r-process sources and give important 
information on the r-patterns produced by these sources. 

Except for the trivial case where an element has a single r-process isotope and some special cases 
|12Uj - jl22j where the hyperfine structure of the atomic spectra allows the extraction of the isotopic 
composition (see discussion of Ba in Sec. I4.3.2|) . stellar observations can only give the total abundance 
of all the isotopes of an element. The observed abundance of an element E in a star is usually given in 
terms of 



where (E/H) is the abundance ratio of the element E relative to hydrogen in the star. Hydrogen is an 
excellent reference element as its overall abundance has changed little over the history of the universe. 
The "metallicity" of a star is measured by 



where the subscript "0" denotes quantities for the sun. There are two possible sources for Fe: SNe II 
and la. However, as SNe la are associated with low-mass progenitors that evolve on timescales of ~ 10^ 
yr, only SNe II with short-lived massive progenitors contributed Fe during the first ~ 10^ yr subsequent 
to the onset of Galactic chemical evolution [t = 0). Over a period of ~ 10^° yr prior to the formation 
of the solar system, SNe II contributed ~ 1/3 of the solar Fe abundance (e.g., ||123ij ). This gives 



for t < 10^ yr. Thus, for stars with [Fe/H] < — 1.5 corresponding to t< 10^ yr, their Fe abundances are 
dominated by SN II contributions. Stars with —3 < [Fe/H] < —2.5 are referred to as ultra-metal-poor 
(UMP) stars and those with —2.5 < [Fe/H] < — 1.5 as metal-poor (MP) stars. 

The observational data of interest here concern the elements Sr and above. In general, both the r- 
process and the main s-process can produce these elements. However, the main s-process contributions 
dominantly come from low-mass AGB stars that evolve on timescales of ~ 10^ yr. So the abundances 
of the elements Sr and above in UMP and MP stars represent the r-process contributions as pointed 
out by Truran |124j . The occurrence of r-process events in the early Galaxy {t < 10^ yr) suggests that 
these events are associated with objects such as SNe II, which have short-lived massive progenitors. 

4.2.1 nonsolar r-patterns in UMP stars 

The observed abundances of a large number of elements in the UMP star CS 22892-052 |1101 1114j are 
shown in Fig. 7. The solar r-pattern [10^ translated to pass through the Eu data is also shown for 
comparison. It can be seen that the data on the heavy r-process elements from Ba to Ir are in excellent 
agreement with the translated solar r-pattern. However, there are large differences between the data 
and this pattern in the region of the light r-process elements below Ba. It is important to examine how 
the uncertainties in the solar r-process abundances may affect the above comparison. The fraction of 
the solar abundance contributed by the s-process is very small {Pq^s ~ 0.3) for the following groups of 



elements: (1) Ru, Rh, and Ag; (2) from Sm to Yb; and (3) Os and Ir (TU] (see Table 1 in Sec. I4.2.2|) . 



loge(E) =log(E/H) + 12, 



(38) 



[Fe/H] = loge(Fe) - loge0(Fe) = log(Fe/H) - log(Fe/H) 



(39) 




(40) 
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Figure 7: The observed abundances in CS 22892-052 (filled circles with error bars: |limill4] ) 
compared with the solar r-pattern (solid curve: JU]) that is translated to pass through the 
Eu data. The data on the heavy r-process elements from Ba to Ir are in excellent agreement 
with the translated solar r-pattern. However, the data on the light r-process elements Rh 
and Ag clearly fall below this pattern. See text for details. See also Sec. 14.2.^ for discussion 
of corrections to the solar r-process abundances of Sr, Y, Zr, Nb, and Ba. 



The solar r-process abundances of these elements should be rather reliable. Thus, in comparing the 
data for CS 22892-052 and the translated solar r-pattern, the agreement for the elements above Ba in 
groups (2) and (3) and the disagreement for Rh and Ag below Ba in group (1) are robust. Regardless of 
the possible uncertainties in the solar r-process abundances of the other elements (see discussion of Sr, 
Y, Zr, Nb, and Ba in Sec. I4.2.2j) . the observed r-pattern in CS 22892-052 is clearly different from the 
overall solar r-pattern. Similar difference is also observed for another UMP star CS 31082-001 |117j . 
These observations provide strong support for the diversity of r-process sources as concluded from the 
meteoritic data on ^'^^Hf and ^^^I. The observed deficiency especially at the light r-process elements 
Rh and Ag relative to the solar r-pattern that is translated to pass through the Eu data reflects the 
characteristics of the H events that produce mainly the heavy r-process nuclei. By inference, a mixture 
of the H events and the L events that produce dominantly the light r-process nuclei is then required 
to account for the overall solar r-pattern. 

4.2.2 yields of H and L events and three-component model for abundances in UMP and 
MP stars 

The differences in the observed r-process abundances between the UMP stars HD 115444 |113j . CS 
31082-001 [mi, and CS 22892-052 [11011111 are shown as Aloge(Z) = logeHDii5444(^)-logecs22892(^) 
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Figure 8: Differences in the observed r-process abundances between HD 115444 jll3j and 
CS 22892-052 (nUl HH (open squares with error bars) and those between CS 31082-001 
(117] and CS 22892-052 (filled circles with error bars). The data are consistent with a 
uniform shift in loge(Z) of ~ —0.6 dex (lower dashed line) for HD 115444 and of ~ 0.21 
dex (upper dashed line) for CS 31082-001, both relative to CS 22892-052. This suggests an 
approximately fixed yield pattern for the H events. The upward offset from the lower dashed 
line for the elements Sr, Y, and Zr in HD 115444 indicates that this star may have received 
large contributions to these elements from the P inventory [to be discussed in the text after 
Eq. (|42j)] in addition to the H events. By contrast, the abundances of Sr, Y, and Zr in CS 
31082-001 and CS 22892-052 are dominated by the contributions from the H events. 



or Iogecs3io82(^) — logecs22892(^) in Fig. 8. It can be seen that within the observational uncertainties, 
these differences represent a uniform shift in loge(Z) of ~ —0.6 dex for HD 115444 and of ~ 0.21 dex for 
CS 31082-001, both relative to CS 22892-052. This means that these three UMP stars have essentially 
the same r-pattern. This pattern is deficient in the light r-process elements such as Rh and Ag relative 
to the solar r-pattern translated to pass through the Eu data (an established result for CS 22892-052, a 
confirmed prediction for CS 31082-001 jll7l 112^, and still a prediction for HD 115444), and therefore, 
is characteristic of the H events (see Sec. I4.2.1|) . The H events then have an approximately fixed yield 
pattern, which can be taken from the data for e.g., CS 22892-052. As the H yield pattern coincides 
with the solar r-pattern and the r-patterns in a number of MP stars |112t I116j in the region of the 
heavy r-process elements, it is reasonable to assume that these elements are produced exclusively by 
the H events. In this case, the occurrence of these events can be represented by the abundance of a 
typical heavy r-process element such as Eu. 

Over the period of Tup ~ 10^° yr prior to the formation of the solar system, a total number 
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= JhTup ~ 10'^ of H events contributed to the solar r-process abundance of Eu for fn ~ (10^ yr) ^. 
This gives (Eu/H)©^^ = n'f^(Eu/R)H, or 

logeH(Eu) = loge0,,(Eu) - logn® ~ -2.5, (41) 

where logeo^^lEu) = 0.52 corresponds to the solar r-process abundance of Eu (TU] and logeH(Eu) 
represents the Eu abundance resulting from a single H event. The data on Eu over the wide range of 
—3 < [Fe/H] < — 1 are shown in Fig. 9. The lowest observed Eu abundances are consistent with the 
enrichment level of logeH(Eu) ~ —2.5 for a single H event. This consistency further establishes that 
the frequency of the H events is fn ~ (10^ yr)~^, which is close to the frequency for replenishing newly- 
synthesized material in the ISM by core-collapse supernovae (see Sec. 14. If fn were ~ (10^° yr)~^, 
which is close to the frequency for replenishing newly-synthesized material in the ISM by neutron star 
mergers (see Sec. 14. HI . n|} would be ~ 1. The corresponding enrichment level of loge/f(Eu) ~ 0.52 for 
a single H event [see Eq. Mlj) ] would be in clear conflict with the Eu data shown in Fig. 9. This again 
suggests that the H events cannot be associated with neutron star mergers |108j . A similar argument 
against the association of the L events with neutron star mergers |108j can be made by using the data 
on the light r-process element Ag ^126] . It is important that these arguments be verified by more 
sophisticated studies on the inhomogeneous chemical evolution of the early Galaxy. 

Figure 9 shows that there is a large dispersion of ~ 2 dex in loge(Eu) over the narrow range of 
-3 < [Fe/H] < -2.5. In particular, CS 31082-001 and CS 22892-052 with the highest Eu abundances 
among the UMP stars shown in this figure have very low values of [Fe/H] = —2.9 and —3.1, respectively. 
This suggests that the H events responsible for Eu enrichment cannot produce any significant amount 
of Fe jl27j . No Fe is expected to be produced in neutron star mergers. However, as argued above, 
the rarity of these events appears to result in conflicts with both the meteoritic data on ^^^Hf and 
observations of abundances in metal-poor stars if they are associated with the H events. While the 
possibility of neutron star mergers being an r-process source merits further studies, the discussion below 
assumes that the H and L events are associated with core-collapse supernovae, which include SNe II 
from the collapse of 0-Ne-Mg and Fe cores as well as the silent supernovae from accretion-induced 
collapse of white dwarfs. A number of SNe II including SN 1987A are known to produce Fe (see e.g.. 
Table 1 of |128j ). In fact, ~ 1/3 of the solar Fe abundance is attributed to SNe II (e.g., |123j ). As no 
Fe can be produced by the H events, this Fe must be assigned to a total number nf = /lTup ~ 10^ of 
L events that occurred with a frequency of /l ~ (10* yr)~^ over the period of Tup ~ 10^° yr prior to 
the formation of the solar system. This gives (l/3)(Fe/H)Q ~ nf (Fe/H)^, or 

[Fe/H]^~-log(3nf)~-2.5, (42) 

where [Fe/H]/, represents the Fe abundance resulting from a single L event. 

The UMP stars have —3 < [Fe/H] < —2.5, and therefore, cannot have received any contributions 
from the L events. The Fe in these stars is attributed to the very massive ( > 100 M©) stars that 
dominated the chemical evolution of the universe until [Fe/H] ~ —3 was reached |ll27j . At this critical 
metallicity, very massive stars could not be formed any more and major formation of regular (~ 1- 
60 Mq) stars took over. In addition to Fe, other elements were also produced by the very massive stars. 
The inventory of such elements associated with the production of [Fe/H] ~ —3 is referred to as the 
prompt (P) inventory. The abundance of an element E in a UMP star represents the sum of the P 
inventory and the contributions from the H events: 

(E/H) = (E/H)p + (E/Eu)^(Eu/H), (43) 

where (E/H)p is the P inventory of E and {E/Eu)h is the relative H yield of E to Eu. The parameter 
(E/Eu)h is assumed to be the same for all H events (see Fig. 8). Then the parameters (E/H)p and 



28 



- 



-0.5 - 



I ' ' ' ' I 

[Fe/H], 



tJO -1.5 

o 



-1 - 



C] • 



-2 - 



-2.5 «3- 



• • 
•••••• • 



log e^(Eu) 



-3 



-2.5 



-2 

[Fe/H] 



-1.5 



-1 



Figure 9: Data on Eu over the wide range of —3 < [Fe/H] < — 1 (filled circles: |115j : open 
squares: |lU9j : asterisks: |113j : filled diamond: |114j : filled triangle: |117p . There is a 
large dispersion of ~ 2 dex in loge(Eu) at —3 < [Fe/H] < —2.5, which suggests that the 
H events responsible for Eu enrichment cannot produce any significant amount of Fe. The 
Eu abundance resulting from a single H event (horizontal dashed line) is consistent with 
the lowest observed Eu abundances. The vertical dashed line indicates the Fe abundance 
resulting from a single L event. See text for details. 



(E/Eu)iif can be obtained from the data for two UMP stars with essentially the same [Fe/H] but very 
different Eu abundances (e.g., CS 22892-052 and HD 115444). It was found that there is a significant 
P inventory for the light r-process elements Sr, Y, and Zr |129j . The parameters (E/H)p and (E/Eu)h 
for these elements are given in terms of logep(E) and logeH(E) in Table 1. The parameter (E/Eu)h 
can be calculated as 

= -|^glogeH(E)-logeH(Eu)^ ^^^^ 

where loge/f(E) represents the abundance of E resulting from a single H event (the H yield of E) and 
logeiy(Eu) = —2.48 corresponds to assigning the solar r-process abundance of Eu to = lO'^ H events 
[see Eq. fl41|) ]. By assuming that the P inventory is negligible for the elements above Zr, the loge/^ 
values for the elements from Nb to Cd can be calculated from the data for CS 22892-052. As the heavy 
r-process elements Ba and above follow the solar r-pattern, the logen values for these elements can be 
calculated from the data for CS 22892-052 or from the solar r-process abundances as in the case of Eu 
[see Eq. Mlj) ]. The logen values for the elements above Zr are also given in Table 1. By applying Eq. 
fl43|) to the data for UMP stars, it was further found that the H events cannot produce any significant 
amount of the elements from Na to Ni, which include not only the Fe group but also the so-called 
"a-elements" such as Mg, Si, Ca, and Ti |130j . 
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Table 1: Parameters of the three-component model for abundances in UMP and MP stars 
(see Table 3 of [129 ). The P inventory in col. (2) and the H yields in col. (3) are calculated 
from the data on HD 115444 HIS] and CS 22892-052 [HI for Sr, Y, and Zr. The P inventory 
for the elements above Zr is assumed to be negligible and set to — oo. The H yields of Nb to 
Cd are calculated from the data for CS 22892-052 and those of Ba to Au from the standard 
solar r-process abundances [10] in col. (5) [corresponding to the solar s-process fractions in 
col. (6)]. The L yields in col. (4) are calculated from cols. (2), (3), and (5) for Sr to Cd 
and set to — oo for Ba and above. The corrected values for Sr, Y, Zr, and Ba are given in 
cols. (7)-(9). These columns are omitted for La to Au, for which elements no corrections 
are recommended. 



E 
(1) 


log ep(E) 
(2) 


loge/^(E) 
(3) 


logei(E) 
(4) 


loge0,^(E) 
(5) 


/?0,.(E) 
(6) 


logen(E) 
(7) 


loge--(E) 
(8) 


/5™7(E) 
(9) 


Fe 


4.5P 


— oo 


5.03^ 












Sr 


0.13 


-1.30 


-oo^ 


1.64 


0.95 


0.35 


2.44 


0.69 


Y 


-1.05 


-2.05 


-1.36 


1.12 


0.92 


-0.42 


1.67 


0.72 


Zr 


-0.13 


-1.53 


-0.40 


1.85 


0.83 


0.26 


2.32 


0.49 


Nb 


— oo 


-2.56 


-2.06 


0.56 


0.85 


c 


c 


c 


Ru 


— oo 


-1.56 


-0.92 


1.60 


0.32 








Rh 


— oo 


-2.36 


-1.20 


1.03 


0.14 








Pd 


— oo 


-1.82 


-0.95 


1.42 


0.46 








Ag 


— oo 


-2.36 


-1.02 


1.14 


0.20 








Cd 


— oo 


-1.91 


-0.85 


1.42 


0.52 








Ba 


— oo 


-1.52^ 


— oo 


1.48 


0.81 


-0.47 


1.79 


0.62 



^ With loge0(Fe) = 7.51 4j, these values correspond to [Fe/H]p = -3 and [Fe/H]i = -2.48. 
^ The standard solar r-process abundance of Sr is saturated by the contributions from 
= 10^ H events, and therefore, does not allow any contributions to Sr from the L events. 
^ The data for BD +17°3248 [118 suggest logei°"(Nb) ^ -0.88 to -0.66, which corresponds 
to loge™;'^(Nb) ^ 1.20 to 1.39 and /^©^.(Nb) ^ 0.37 to 0.02. This corrected L yield of Nb 
should be tested by future data for MP stars. 

This H yield of Ba, which is calculated by attributing its standard solar r-process abun- 
dance to n% = 10^ H events, essentially coincides with the value of loge//(Ba) ~ —1.57 
obtained from the data for CS 22892-052. The use of the latter value is recommended. 



30 



Table 1 (continued) 



E 


log ep(E) 


loge/^(E) 


logeL(E) 


log 60,,. (E) 


/50,.(E) 








[V 






La 


— oo 


-2.22 


— oo 


0.78 


0.62 


Ce 


— oo 


-2.02 


— oo 


0.98 


0.77 


Pr 


— oo 


-2.51 


— oo 


0.49 


0.49 


Nd 


— oo 


-1.90 


— oo 


1.10 


0.56 


Sm 


— oo 


-2.20 


— oo 


0.80 


0.29 


Eu 


— oo 


-2.48 


— oo 


0.52 


0.06 


Gd 


— oo 


-2.00 


— oo 


1.00 


0.15 


Tb 


— oo 


-2.70 


— oo 


0.30 


0.07 


Dy 


— oo 


-1.92 


— oo 


1.08 


0.15 


Ho 


— oo 


-2.53 


— oo 


0.47 


0.08 


Er 


— oo 


-2.13 


— oo 


0.87 


0.17 


Tm 


— oo 


-2.93 


— oo 


0.07 


0.13 


Yb 


— oo 


-2.22 


— oo 


0.78 


0.33 


Lu 


— oo 


-2.98 


— oo 


0.02 


0.20 


Hf 


— oo 


-2.61 


— oo 


0.39 


0.56 


A-Ci 




o. ou 




—0 "^fi 


41 


w 


— oo 


-2.68 


— oo 


0.32 


0.56 


Re 


— oo 


-2.75 


— oo 


0.25 


0.09 


Os 


— oo 


-1.66 


— oo 


1.34 


0.09 


Ir 


— oo 


-1.63 


— oo 


1.37 


0.01 


Pt 


— oo 


-1.34 


— oo 


1.66 


0.05 


Au 


— oo 


-2.20 


— oo 


0.80 


0.06 
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The Fe at —2.5 < [Fe/H] < — 1.5 represents the contributions from the L events added to the P 
inventory. The abundance of an element E in an MP star in this metalhcity range can be written as 



(E/H) = (E/H)p + (E/Eu)^(Eu/H) + (E/Fe)^[(Fe/H) - (Fe/H)p], 



(45) 



where (E/Fe)z, is the relative L yield of E to Fe and is assumed to be fixed for all L events. With the 
parameters (E/H)p and (E/Eu)^ given in terms of logep(E) and loge//(E) in Table 1, the parameter 
(E/Fe)^ can be obtained from the solar r-process abundance of E: 



where 1/3 of the solar Fe abundance is taken to be contributed by the L events. The parameter (E/Fe)i 
calculated in this way for the light r-process elements from Sr to Cd are given in terms of logeL(E) in 
Table 1. The parameter (E/Fe)^ can be calculated as 



where loge/,(E) represents the abundance of E resulting from a single L event (the L yield of E) and 
logeL(Fe) = 5.03 ([Fe/Hj^ = —2.48) corresponds to assigning 1/3 of the solar Fe abundance to = 10^ 
L events [see Eq. (jH))]. With the parameters logep(E), logeff(E), and logeL(E) given in Table 1, the 
abundance of E in any UMP or MP star can be calculated from Eqs. ()43p - ()45|) and ()47|) by using only 
the observed abundances of Eu and Fe in the star. This is the three-component (P, H, and L) model 
for abundances in UMP and MP stars as developed by Qian and Wasserburg |129j . 

By comparing the data and the model predictions, large discrepancies were found for the elements 
Sr, Y, Zr, and Ba at —2.5 < [Fe/H] < — 1.5. One possible explanation is that significant s-process 
contributions to these elements already occurred at metallicities as low as [Fe/H] ~ —2.5 (e.g., |115j ) 
but were not taken into account by the three-component model. However, it is also possible that the 
contributions from the L events to the elements Sr, Y, and Zr were greatly underestimated by the model 
and that the L events may produce a large amount of Ba. Good agreement between the data and the 
model can be obtained when the solar r-process abundances of Sr, Y, Zr, and Ba are increased from their 
standard values by factors of ~ 2-6 and the L yields of these elements are correspondingly increased 
|129j . The solar abundances of Sr, Y, Zr, and Ba are dominated by the s-process contributions (/?0^s>O.8 
[TUj : see Table 1). Thus, the solar r-process abundances of these elements obtained by subtracting the 
s-process contributions are subject to possibly large uncertainties (see Sec. II. 1|) . It remains to be seen if 
future s-process studies can justify the corrected solar r-process abundances obtained by comparing the 
data and the three-component model for abundances in UMP and MP stars. The corrected parameters 
are given in Table 1. The corrected L yield of Ba is ~ 10 times larger than the H yield of Ba. However, 
the heavy r-process elements above Ba are still produced exclusively by the H events. This suggests 
that the L events can only produce the elements up to and including Ba. 

As an example to test the three-component model (with corrected L yields of Sr, Y, Zr, and Ba), 
the abundances calculated for BD -|-17°3248 by using the observed values of loge(Eu) = —0.67 and 
loge(Fe) = 5.43 are compared with the data |118j in Fig. 10. It can be seen that the agreement 
between the model and the data is quite good except for the element Nb. However, as in the case of 
the elements Sr, Y, Zr, and Ba, the solar abundance of Nb is dominated by the s-process contributions 
{(3q^s = 0.85 lOj see Table 1). Thus, the L yield of Nb derived from Eq. (j46|) by using its solar 
r-process abundance may be quite uncertain. The data for BD -|-17°3248 suggest a corrected L yield of 
loge™'''^(Nb) ^ —0.88 to —0.66, which corresponds to a rather small value of /50°"(Nb) ^ 0.37 to 0.02. 
These corrected values should be tested by future data for MP stars and detailed s-process models. 



(E/H)^,, = (E/H)p + (E/Eu)^(Eu/H)q^ + (E/Fe) Jl/3)(Fe/H) 



0' 



(46) 




(47) 
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Figure 10: The observed abundances in BD +17°3248 (filled circles with error bars: |118j ) 
compared with the predictions from the three-component model (solid curves). The model 
|129j uses only the data on Eu and Fe to calculate the abundances of all the other elements. 
The agreement between the model and the data is quite good except for the element Nb. 
The L yield of Nb may have been underestimated by the model due to the uncertainty in 
the solar r-process abundance of Nb. See text for details. 



4.2.3 variations in the H yield pattern and nucleochronology 

The discussion of stellar observations has focused on the stable elements so far. Observations of the 
long-lived radioactive elements Th and U in metal-poor stars are considered below in connection with 
nucleochronology. The dominant radioactive nuclei in stars with an age of tgtar ~ 10^'' yr are ^^^Th (with 
a lifetime of r232 = 2.03 x 10^° yr) and ^ss^j [j^^^ = g.45 x 10^ yr). The period over which the r-process 
events could contribute to the initial abundances of these two nuclei in such old stars is limited by 
the age of the universe (?« 1.5 x lO^'' yr) to shorter than the lifetime of either nucleus. So the initial 
abundance ratio (^^®U/^^^Th)o in these stars essentially reflects the relative yield of ^38^ 232rj]^ ^^^^ 
the r-process events. As these two nuclei are close in mass number, their relative yield V238/V232 may 
be rather constant among individual events. Then the present abundance ratio (^^^U/^^^Th)star in an 
old star is 

which gives 

W ~ (21.8 X 10^ yr) [log(F238/>^232) - log(U/Th),tar] , (49) 
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where the present abundance ratio (^'^^U/^'^^Th)star in the star has been set to be the same as the 
elemental abundance ratio (U/Th)star- In principle, 1238/^232 can be calculated from theory and various 
estimates give ^238/^232 ~ 0.53 ± 0.23 (e.g., |131j ). If the abundances of Th and U are measured in 
the star, then its age can be estimated from Eq. So far, the abundances of both Th and U 

have been observed only in one UMP star, CS 31082-001 |II71[IH21, and one MP star, BD +17°3248 
jll8j . Age estimates based on Eq. ()49|) are rather unreliable for BD +17°3248 mainly due to the large 
observational uncertainty of ~ 0.3 dex in loge(U) for this star. By contrast, an age of (13.5 ±5.0) x 10^ 
yr can be estimated from log(U/Th) = -0.94 ± 0.11 for CS 31082-001 [mj (see also [02]), where the 
error in the age is dominated by the uncertainty in the theoretical estimate of ^233/ ^232- 

If the heavy r-process elements including the long-lived chronometers are exclusively produced by 
the H events and the H yield pattern is fixed, then nucleochronology can be applied to stars for which 
the abundances of Th and a stable heavy r-process element such as Eu are known (e.g., |113|, lllbl I118| 
11331 1134j ). The present abundance ratio (Th/Eu)star in a star is related to its age by 

/Th\ / tstar 

star 




which gives 

tstar = (46.7 X 10^ yr)[log(Th/Eu)o - log(Th/Eu)star]. (51) 

The initial abundance ratio (Th/Eu)o in Eqs. fl50|) and (j^T|) is essentially the relative H yield of Th 
to Eu, (Th/Eu)i^, and can be estimated from the r-process abundance ratio (Th/Eu)Ess,r in the early 
solar system. Due to the long lifetime of ^^^Th, the replenishment of this nucleus in the ISM by the H 
events can be considered as continuous and the change in (Th/H) is governed by 

d /Th\ , /Th\ 1 /Th\ , , 



dt\n •' V H 



H 



which can be solved to give 




ESS 




exp 



Tup \ 

^232 / 



(53) 



where Tup is the period of uniform production prior to the formation of the solar system. The abundance 
ratio (Th/Eu)Ess,r is then 




ESS,r 




1 — exp 



T, 



UP 



T232 



(54) 



where (Eu/H)Ess,r = (Eu/H)^/hTup has been used. With (Th/Eu)Ess,r = 0.042/0.0917 111101, Eq. 

gives (Th/Eu)ff = 0.58 for Tup = 10^° yr. The value of (Th/Eu)/^ only decreases slightly to 0.54 
if Tup = 7 x 10^ yr is used instead. 

For CS 22892-052, log(Th/Eu) = -0.67 ±0.11 Applying Eq. ^ with (Th/Eu)o = 

(Th/Eu)^ = 0.58 gives an age of (20.1 ± 5.1) x 10^ yr for this star (see also [1331 CMl). However, the 
same procedure applied to CS 31082-001 with log(Th/Eu) = — 0.22±0.05 J117| leads to an unacceptable 
age of (-0.9 ±2.3) x 10^ yr. Clearly, Eq. ^ must be applied to CS 22892-052 and CS 31082-001 with 
different values of (Th/Eu)o. By using the age of (13.5 ±5.0) x 10^ yr for CS 31082-001 estimated above 
from the observed abundances of Th and U in this star, a value of (Th/Eu)o = 1.20 ± 0.26 is obtained 
from Eq. (j3T|) . The value of (Th/Eu)/^ = 0.58 inferred above from Eq. ([3^ should only represent some 
average over many H events. The data for CS 31082-001 indicate that the relative yield of Th to Eu 
in individual H events can vary by a factor of ~ 2 around this average. Such variations not only have 
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severe impact on the Th/Eu chronology (e.g., |135j ). but also have profound implications for r-process 
nucleosynthesis (e.g., |l36j). Fission cycling tends to result in a global steady flow between the fission 
fragments and the fissioning nuclei, and hence, a robust yield pattern of the heavy r-process elements 
Ba and above (e.g., [SI])- This cannot occur during all the H events in view of the variations in the 
relative yield of Th to Eu. The possible role of fission (but not fission cycling) in producing a generic 
H yield pattern and the associated variations is discussed in Sec. 14.3.21 

4.3 Observational implications for the r-process 

The meteoritic data on ^^^Hf and ^^^I require diverse sources for the r-process and are consistent with 
the replenishment of these two nuclei in the ISM by the H and L events at frequencies of ~ (10^ yr)~^ 
and ~ (10^ yr)~^, respectively. These frequencies may be explained by associating both kinds of events 
with core-collapse supernovae (see Sec. 14.11) . The solar r-process abundances resulted from ~ 10^ H 
events and ~ 10^ L events. The corresponding enrichment level by a single H or L event is in agreement 
with the observed r-process abundances in metal-poor stars (see Fig. 9 and 108J). The mass fraction 
of the heavy or light r-process elements in the solar system is X^^^^^ ~ ^©^r* ~ 4 x 10~^ JU]- If the 
amount of ISM that will mix with the ejecta from each event is Mmix ~ 3 x 10^ Mq as in the case 
of supernovae, XQ^^^-^Mmix ~ -^©^r^^mix ~ 10~^ Mq of the heavy or light r-process elements must be 
provided by the ~ 10'^ H events or the ~ 10^ L events that occurred prior to the formation of the solar 
system. This requires ~ 10~^ Mq or ~ 10~^ M© of the r-process material be ejected by each H ot L 
event, respectively. The neutrino-driven wind from a protoneutron star can eject ~ 10~^-10~^Mo of 
material, although it remains to be established that the conditions required for the r-process can be 
obtained in the wind (see Sec. 13.11) . The mass fraction of Fe in the solar system is 10~^ j3]. In 

order to account for ~ 1/3 of this Fe by the ~ 10^ L events that occurred prior to the formation of the 
solar system, each L event must eject ~ {X^ /300)Mjnix ~ 0.1 Mq of Fe. This is in agreement with the 
amount of Fe inferred from the light curves for a number of SNe II (see e.g.. Table 1 of |128j ). The rarity 
of neutron star mergers appears to result in conflicts with both the meteoritic data and observations of 
abundances in metal-poor stars if they are associated with the H or L events. The arguments leading 
to this conclusion must be examined by more sophisticated models of Galactic chemical evolution that 
include the details for the mixing of the ejecta from individual nucleosynthetic events with the ISM. 
While the possibility of neutron star mergers being an r-process source should also be studied, the 
discussion below again assumes that the H and L events are associated with core-collapse supernovae. 

4.3.1 three-component model for abundances in UMP and MP stars and evolution of Eu 
abundance relative to Fe 

In the three-component model for abundances in UMP and MP stars (see Sec. I4.2.2j) . the enrichment 
level for a single H ot L event corresponds to loge//(Eu) ~ —2.5 or [Fe/H]^ ~ —2.5. For a mixture of 
a number uh of H events and a number hl oi L events in the ISM, the corresponding abundances of 
Eu and Fe are given by 



where [Fe/H]p ~ —3 represents the P inventory of Fe. On average, a fraction q = fn/ifH + /l) of the 
r-process events is of the H kind (g ~ 0.9 for the three-component model). The probability for an ISM 
to have a mixture of a number hh of H events and a number hl oi L events is 




(55) 
(56) 



P{nH,nL) 



{nH + nL)\ 



g"«(l-g) 



riL+l 



(57) 
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Figure 11: Data on Eu at —3 < [Fe/H] < — 1 (same symbols as in Fig. 9) compared with the 
99% confidence contours (sohd curves) and the mean trend (dashed curve) calculated from 
the three-component model |137j . There is general agreement between the model and the 
data. The several data below the lower 99% confidence contour at [Fe/H] > —1.5 can be 
explained by taking into account the Fe contributions from SNe la. The extremely high Eu 
abundances in CS 31082-001 and CS 22892-052 are very close to the upper 99% confidence 
contour and may not refiect the chemical composition of the ISM from which these two stars 
formed. Instead, such high r-process enrichments may be due to surface contamination by 
the ejecta from core-collapse supernovae {H events) in binaries |125j . See text for details. 



The 99% confidence contours for log e(Eu) over a range of [Fe/H] and the mean trend for evolution of Eu 
abundance relative to Fe can be calculated from Eqs. The results for loge//(Eu) = —2.48, 

[Fe/H]p = —3, [Fe/H]j^ = —2.48, and q = 10/11 |137j are compared with the data for metal-poor 
stars in Fig. 11. It can be seen that essentially all the data lie within the calculated 99% confidence 
contours. The data are also in general agreement with the calculated mean trend. There are several 
data below the calculated lower 99% confidence contour at [Fe/H] > —1.5 (see Fig. 11). This is because 
the calculations assumed that only the L events added Fe beyond the P inventory and did not take into 
account the SN la contributions to Fe at [Fe/H] > —1.5. Therefore, the observed loge(Eu) values over 
the wide range of —3 < [Fe/H] < — 1.5 are in accord with the three-component model for abundances in 
UMP and MP stars. Note that essentially all the other models (e.g., |138| I139j ) focused on explaining 
[Eu/Fe] = log(Eu/Fe) — log(Eu/Fe)0, but not loge(Eu), as a function of [Fe/H]. 

The data for the UMP stars CS 31082-001 and CS 22892-052 lie very close to the calculated 
upper 99% confidence contour for loge(Eu) in the region where only the H events could occur. The 
observed value of loge(Eu) = -0.76 in CS 31082-001 Hm or loge(Eu) = -0.93 in CS 22892-052 mH 
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corresponds to uh = 52 or 35 H events [logeiif(Eu) = —2.48]. With a frequency ratio of fu/ fi = 10, 
it is unlikely that these many H events could have occurred in an ISM without any accompanying L 
events. Thus, the observed r-process abundances in CS 31082-001 and CS 22892-052 may not reflect the 
composition of the ISM from which they formed. Instead, each star may have had a binary companion, 
which evolved and exploded as a supernova H event. The ejecta from the explosion then contaminated 
the surface of CS 31082-001 or CS 22892-052 with high r-process enrichments |125j . This scenario has 
two possible outcomes. If the binary survived the explosion, then the highly r-process enriched star 
would have a neutron star or black hole companion and exhibit periodic shifts in its radial velocity. 
Such shifts were indicated by observations of CS 22892-052 |14()j . On the other hand, if the binary was 
disrupted by the explosion, then the highly r-process enriched star might acquire a large space velocity 
from its previous orbital motion. This may be revealed by proper motion measurements of the star. 

4.3.2 fission and the H yield pattern 

It is important that observations be carried out to test the above surface contamination scenario for 
explaining the highly r-process enriched metal-poor stars. This will not only help establish core-collapse 
supernovae as the r-process site, but also help determine the yield patterns of individual supernovae. If 
both CS 31082-001 and CS 22892-052 are cases of surface contamination, then the observed r-patterns 
in these two stars represent the yield patterns of two individual H events. Figure 8 shows that these 
two patterns are essentially the same in the region of the elements from Sr to Ir. However, the relative 
H yields of Th to Eu inferred from the data for CS 31082-001 and CS 22892-052 are very different (see 
Sec. I4.2.3|) . Thus, the H yield pattern has some generic feature but also exhibits certain variations. 
Concerning the generic H yield pattern, there is a lack of data between Cd {A = 111-114, 116) and 
Ba {A = 135, 137, 138). However, the meteoritic data require that the H events produce very little 
^^^I (see Sec. 14. 1|) . Observations also indicate that the Ba in the UMP star HD 140283 is dominantly 
^^^Ba. This determination of the Ba isotopic composition is possible because only the odd-A isotopes 
^^^Ba and ^'^''Ba have finite nuclear magnetic moments that give rise to hyperfine structure of the Ba 
spectra. The fraction of the odd-A Ba isotopes in HD 140283 was found to be /odd = 0.08 ± 0.06 |120j 
and 0.30 ± 0.21 [122, . Thus, there appears to be a rather sharp rise in the H yield at ^'^^Ba. This is 
consistent with the very low H yields at A ~ 130 required by the meteoritic data. 

The above discussion of the H yield pattern raises two issues: how can the light {A < 130) and 
heavy [A > 130) r-process elements be produced without any significant production at A ~ 130 and 
how can some variations in the yield pattern be accommodated? Both issues may be resolved by 
considering fission of the progenitor nuclei during decay towards stability after the r-process freezes out 
in core-collapse supernovae |136j . It is assumed that fission does not occur during the r-process (i.e., 
no fission cycling) and that the r-process produces a freeze-out pattern covering 190 < A < 320 with 
a peak at A ~ 195. As the progenitor nuclei decay towards stability after the freeze-out, all of those 
with 260 < A < 320 eventually undergo spontaneous fission (e.g., |Hn|)- Due to the strong influence of 
the closed proton and neutron shells at ^^^Sn, the fission of the progenitor nuclei with 260 < A < 320 is 
expected to produce one fragment at A ~ 132 and the other at 130 < A < 190. Some of the progenitor 
nuclei with 230 < A < 260 would also undergo spontaneous fission during decay towards stability, thereby 
producing one fragment again at A ~ 132 and the other at 100 ^ A < 130. The possibility of fission 
for the progenitor nuclei with 230 < A < 260 may be greatly enhanced by reactions with the neutrinos 
emitted in core-collapse supernovae as these nuclei could be highly excited by such reactions (e.g., |52]). 
Neutrino reactions may even induce fission of the progenitor nuclei with 190 < A < 230. Experiments 
using energetic particles to induce fission of the stable or long-lived nuclei in this mass range showed 
that the fission mode is dominantly symmetric with no preference for a fission fragment at A ~ 132, 
and the mass ratio of the two fission fragments is ~ 1-1.2 (e.g., |14HI142] ). So neutrino-induced fission 
of the progenitor nuclei with 190 ^ A < 230 is expected to produce fragments at 86 < A < 125. 
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In the above r-process scenario |136j . nuclei with 86< A < 190 would be produced by fission of the 
progenitor nuclei with 190 < A < 320 during decay towards stability after the freeze-out. The fission 
yields may be modified by capture of the neutrons released from fission. For example, the accumulation 
of the fission fragments at A ~ 132 may be shifted to somewhat higher A. This may explain the sharp 
rise in the H yield at ^^^Ba indicated by the isotopic composition of Ba in HD 140283. In any case, there 
would be very little production of the nuclei with A ~ 130. The ratio of the fission yields at 8Q<A < 190 
relative to the surviving abundances at 190 < A < 260 depends on the number of neutrino reactions 
experienced by each progenitor nucleus after the freeze-out. Variations of this number among individual 
supernovae may explain the significant difference in the abundance ratio Th/Eu between CS 31082-001 
and CS 22892-052. The progenitor nuclei far from stability deexcite mainly by neutron emission. The 
probability for deexcitation by fission increases as the charge of the progenitor nuclei increases during 
decay towards stability. Thus, neutrino-induced fission is closely associated with neutrino-induced 
neutron emission. The level of neutrino interaction (a few neutrino reactions per progenitor nucleus) 
after the freeze-out required to account for the H yield patterns is consistent with that for producing 
the nuclei with A = 183-187 by neutrino-induced neutron emission from the progenitor nuclei in the 
yield peak at A ~ 195 (see Fig. 5) |136j . 

4.3.3 supernova progenitors for H and L events 

The three-component model for abundances in UMP and MP stars (see Sec. I4.2.2|) requires that the H 
events produce very little of the elements from Na to Ni including Fe |130j . It is also argued above that 
these events are associated with core-collapse supernovae, which include SNe II from the collapse of 
0-Ne-Mg and Fe cores as well as the silent supernovae from accretion-induced collapse (AIC) of white 
dwarfs. It is commonly thought that in SNe II, the elements below Si are produced by hydrostatic 
burning during presupernova evolution and those from Si to Ni are mostly produced by explosive 
burning associated with the shock propagation. In order not to produce a significant amount of the 
elements from Na to Ni, an H event must either have a presupernova structure lacking substantial 
hydrostatic burning shells or for some reason, have all the material below the He burning shell fall 
back onto the central remnant |130j . Stars with masses of ~ 8-10 Mq have very thin shells at the end 
of their lives and explode due to the collapse of an 0-Ne-Mg core instead of an Fe core These 
stars could be the progenitors for the H events. A supernova from 0-Ne-Mg core collapse would only 
produce ~ 0.002 of Fe jl21- This is ~ 50 times smaller than the Fe yield of ~ 0.1 for an L 
event. In addition, the AIC of a bare white dwarf will not produce any of the elements from Na to Ni, 
and therefore, can also be associated with the H events. A small amount of material would be ejected 
in the neutrino-driven wind from the protoneutron star produced by any core collapse. This wind is a 
possible r-process site although it remains to be seen if adequate conditions for the r-process can be 
obtained in the wind (see Sec. 13. ip . Thus, both SNe II from 0-Ne-Mg core collapse and AIC of white 
dwarfs may correspond to the H events or at least a subset. Another possibility to accommodate the 
nucleosynthetic requirements of the H events is that SNe II from Fe core collapse of progenitors at the 
higher-mass end would suffer severe fall back |143j . It has been argued that SN 1997D, which ejected 
only 0.002 M0 of Fe, was such a case |144[ I145j . A potential problem with such a scenario for the H 
events is that the severe fall back may be an obstacle for ejection of the r-process material from the 
inner most regions of the supernova. 

As mentioned earlier, the Fe yield of ~ 0.1 assigned to the L events is consistent with the 
amount of Fe inferred from the light curves for a number of SNe II including SN 1987A (see e.g.. Table 
1 of |128j ). In addition, overabundances of Sr and Ba were observed in the spectra of SN 1987A |14tij - 
|148j . These were attributed to production by the weak s-process during core He burning prior to the 
explosion |149j . However, the observed abundance ratio Ba/Sr is too large for the weak s-process to 
explain. This favors an r-process origin for the observed overabundances of Sr and Ba in SN 1987A 
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|15Uj . By assuming a mass distribution identical to that of ^^Ni inferred from the spectra and the hght 
curve, it was estimated that SN 1987A produced ~ 6 x 10~^ Mq of Ba |15Uj . For a dilution mass of 
~ 3 X 10^ Mq, this Ba yield corresponds to loge(Ba) ~ 0.2 and is of the same order of magnitude as 
the Ba yield assigned to an L event [logeL(Ba) = —0.47; see Table 1]. The observed abundance ratio 
Ba/Sr ~ 2.5(Ba/Sr)0 in SN 1987A [HZI is also consistent with the L yield ratio (Ba/Sr)^ ~ (Ba/Sr)© 
(see Table 1 and The progenitor of this supernova was a 20 M© star (e.g., |151j ). The progenitors 
of SN 1993J [ina CnH] and SN 19941 [1111 with approximately the same Fe yields as that of SN 1987A 
were ~ 12-15 Mq stars. The progenitor of SN 1997D with an extremely low Fe yield of 0.002 Mq might 
be a 26 star |144j . It is possible that the L events are associated with SNe II from progenitors of 
~ 12-25 Mq. The SNe II from progenitors outside this mass range, as well as AIC of white dwarfs, may 
then be associated with the H events. 

4.4 Other r-process studies based on abundances in metal-poor stars 

While only core-collapse supernovae and neutron star mergers are discussed here as the possible r- 
process site, many alternatives were proposed in the past (see jH] for a review). The various candidate 
sites differ in the evolution timescale of the associated astrophysical object (e.g., SNe II have short-lived 
progenitors) or in the dependence on the metallicity of the associated astrophysical environment (e.g., 
some proposed sites rely on the initial metal abundance for the seed nuclei or the neutron source). 
These differences would lead to different evolution of the r-process abundances in the Galaxy. Thus, 
it is possible to discriminate the proposed r-process sites based on observations of stellar abundances 
over a wide range of metallicity |155j . Most studies used data on [E/Fe] as a function of [Fe/H], 
where [E/Fe] = log(E/Fe) — log(E/Fe)0 with E being an element such as Eu. Interpreting the data 
at [Fe/H] < —2.5 as indicating an increase of [Eu/Fe] with increasing [Fe/H], some studies |155| 1156] 
claimed that Eu production was delayed compared with Fe production in the early Galaxy. Assuming Fe 
production by most SNe II, these studies concluded that the r-process occurs in SNe II from longer-lived 
progenitors at the lower- mass end (~ 8-10 M©). With more data (see Fig. 9), it is clear that there is a 
large dispersion in [Eu/Fe] rather than a definite trend at — 3< [Fe/H] < —2.5. This dispersion is due to 
inhomogeneous chemical evolution in the early Galaxy, where the abundances in an ISM were strongly 
influenced by a small number of nucleosynthetic events. Thus, although SNe II from progenitors of 
~ 8-10 Mq may be an r-process site, it is inappropriate to justify this by assuming a smooth trend for 
evolution of [Eu/Fe] at low [Fe/H]. 

The inhomogeneous nature of chemical evolution in the early Galaxy was recognized by many 
authors (e.g., |lll|,ll38lll39j and |157j - [TB^ ). Some authors |,159j followed Galactic dynamics and metal 
enrichment with an iV-body smooth particle hydrodynamic code, in which nucleosynthetic products 
from a stellar particle were transferred to the neighboring gas particles according to the distance between 
the stellar particle and the gas particle. The observed large dispersion in [Ba/Fe] at [Fe/H] < —2.5 
similar to that in [Eu/Fe] was reproduced by this approach. The observed large dispersions in [Sr/Fe], 
[Ba/Fe], and [Eu/Fe] at [Fe/H] < —2.5 were also accounted for by a different approach that treated 
the chemical evolution of the Galactic halo in terms of a series of cloud coalescence and fragmentation 
|161j . Yet another approach focused on the role of individual supernovae. By assuming that stars form 
inside supernova remnants with a chemical composition reflecting the mixture of the supernova ejecta 
and the ISM swept by the remnant, it was shown that the large dispersion in [Eu/Fe] at [Fe/H] < —2.5 
could be explained |138| 11391 1160j . Based on this approach, it was found that the data are consistent 
with assigning Fe production to SNe II from progenitors of > 10 Mq and Eu production to those from 
progenitors of either 8-10 M© or > 30 M© |138j . In any case, not all SNe II can produce both Fe and 

Eu nsi. 

A more extreme approach assumed that the abundances in stars with [Fe/H] < —2.5 represent the 
yields of individual supernovae |158j (see also |157j ). The progenitor of an SN II responsible for the 
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abundances in a star was identified by comparing the observed [Mg/H] with the theoretical SN II yields 
of Mg, which are more or less consistently calculated by different models. The observed abundances 
of the other elements in the star were then used to infer the corresponding yields for the responsible 
SN II. The inferred Fe yields were in agreement with the amount of Fe estimated from the light curves 
for a number of SNe II |158j . In addition, it was found that the SN II yields of Mn, Fe, Co, and Ni 
inferred from the abundances in metal-poor stars increase while that of Eu decreases with increasing 
progenitor mass |163j . The observed large dispersion in [Ba/Fe] at —3 < [Fe/H] < —2.5 was attributed 
to a sharp decrease in the SN II yield of Ba from 8.5 x 10~^ Mq for a 20 Mq progenitor to 4.5 x 10~^ Mq 
for a 25 M0 progenitor |l()4j . Furthermore, the observed overabundances of Sr and Ba in SN 1987A 
(see Sec. I4.3.3|) were used to limit the progenitor masses to 20 ± 0.7 for SNe II producing Ba and 
possibly 20 lb 3 Mq for those producing Sr |150j . While these results are extremely quantitative and the 
underlying approach is rather interesting, some problems with the adopted interpretation of the data for 
metal-poor stars cannot be ignored. For example, the claimed sharp decrease in the SN II yield of Ba 
with increasing progenitor mass was based on interpreting the data on [Ba/Fe] at —3 < [Fe/H] < —2.5 
as indicating a strict decrease of [Ba/Fe] with increasing [Fe/H]. However, in consideration of the 
observational uncertainties in [Fe/H], no clear trend except for a large dispersion can be identified from 
e.g., the observed values of [Ba/Fe] = 0.9, 0.18, and 1.17 for CS 22892-052 [nU], HD 115444 dig, and 
CS 31082-001 mZI with [Fe/H] = -3.1, -2.99, and -2.9, respectively. 

All the above approaches to understand abundances in metal-poor stars merit attention. In ad- 
dressing the observed large dispersions in r-process abundances at —3 < [Fe/H] < —2.5, most of these 
approaches assumed that Fe is produced by the majority of SNe II while the r-process elements are 
produced by a small subset. In contrast, the three-component model discussed in Sec. 14.2.21 attributed 
the Fe at —3 < [Fe/H] < —2.5 to the P inventory that was dominantly produced by the very massive 
stars prior to the onset of major formation of regular stars and assigned the subsequent Fe produc- 
tion by SNe II to the low- frequency L events. The large dispersion in the abundance of e.g., Eu at 
—3 < [Fe/H] < —2.5 was then caused by the occurrence of the high-frequency H events producing Eu 

but no Fe [HZI. Some of the stars with [Fe/H] 3 such as CS 31082-001 and CS 22892-052 were 

argued to be cases of surface contamination by the ejecta from the {H event) explosion of a binary 
companion |125j . Future observations of more stars with extreme r-process enrichments but low [Fe/H] 
along with the tests of the surface contamination scenario discussed in Sec. 14.3.11 should help identify 
the correct approach to understand abundances in metal-poor stars. 

5 Summary and Outlook 

There have been many advances in theoretical and observational studies of the r-process since the last 
major review by Cowan et al. jH] in 1991. This review tries to cover most of these advances. A brief 
summary is given below with some comments on possible further progress. 

The conditions required for the r-process have been examined in terms of the electron fraction, the 
entropy per baryon, and the dynamic timescale for adiabatic expansion from high temperature and 
density (see Sec. EI). This parametrization of the astrophysical environment treats the determination of 
the neutron-to-seed ratio, the r-process, and the freeze-out in a smooth sequence and is more realistic 
than the traditional parametrization using a set of neutron number density, temperature, and neutron 
irradiation time. The questions to be addressed are: how will the conditions required for the r-process 
change if this process occurs in an intense flux of supernova neutrinos (see Sec. 13.1.11) and is it possible to 
refine these conditions by requiring the r-process to produce the yield patterns inferred from meteoritic 
data and observations of metal-poor stars (see Sec. l4.2.!Zj) ? 

There have been many developments regarding core-collapse supernova and neutron star merger 
models of the r-process (see Sec. Ej). In particular, the neutrino-driven wind model associated with 
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core-collapse supernovae has been scrutinized in substantial detail. While current models do not give 
the conditions for an r-process in the wind, meteoritic data on ^^^Hf and ^^^I as well as observations of 
metal-poor stars are in support of core-collapse supernovae being the r-process site (see Sec. EJ. The 
wind model is attractive in that neutrino interaction can eject the right amount of material as required 
from each supernova to account for the solar r-process abundances. The essential role of neutrinos 
in this model also provides a connection between r-process nucleosynthesis and neutrino properties 
such as flavor transformation. The crucial problem is: what remedies (see Sec. I3.1.2|l can lead to an 
r-process in the wind? As for the neutron star merger model, future studies should try to reduce the 
uncertainties in the amount and the neutron-richness of the ejecta. The neutrino-driven wind associated 
with neutron star mergers (see Sec. 13. 3|) is worth investigating in detail. The observational consequences 
of the neutron star merger model, especially the challenges from meteoritic data and observations of 
metal-poor stars (see SecE}, need to be addressed. 

The evidence from meteoritic data and observations of metal-poor stars for the diversity of r- 
process sources (see Sec. E} represents a breakthrough in r-process studies. The demonstration of 
some regularity in r-patterns and of large dispersions in r-process abundances at low [Fe/H] by stellar 
observations has had profound impact on the understanding of r-process nucleosynthesis (see Sec. |3]). 
Future observations should study in more detail the regularities and the variations of r-patterns in 
metal-poor stars. In this regard, it is extremely important to study the abundance ratios of the other 
heavy r-process elements relative to Eu, such as Ba/Eu, La/Eu, and Nd/Eu, at — 2.5< [Fe/H] < — 1.5. If 
the heavy r-process elements are indeed produced by a unique kind of events with a fixed yield pattern, 
then these ratios should stay constant until the onset of the s-process contributions at [Fe/H] ~ —1.5. 
Extensive data are currently available only for Ba, La, Nd, and Dy and exhibit rather constant values 
of La/Eu and Dy/Eu, large spread in Nd/Eu with no clear trend, and a systematic increase of Ba/Eu 
with increasing [Fe/H] (see Fig. 4 of |116j and Fig. 3 of jl29j ). The behavior of Ba/Eu and La/Eu at 
— 2.5<[Fe/H]< — 1.5 indicates large contributions to Ba but no contributions to the next element La from 
another kind of r-process events associated with Fe production (see Sec. I4.2.2|) . This needs to be tested 
by more extensive and precise data in the future. More data on Th and U in metal-poor stars would 
be of great value to studies of variations in the yield pattern of the heavy r-process elements (see Sec. 
I4.2.3|) . A substantial data base for abundance patterns in very metal-poor stars with extreme r-process 
enrichments may provide unique opportunities to study the r-process. Measurements of radial velocity 
and proper motion for such stars are crucial in testing the scenario of surface contamination by the 
ejecta from core-collapse supernovae in binaries (see Sec. 14. 3. This is the most promising approach to 
identify core-collapse supernovae as the r-process site, although a more direct approach is to detect the 
7 rays from the decay of r-process progenitor nuclei in supernova remnants (e.g., |165| 116b'] ). Another 
approach to demonstrate the association of core-collapse supernovae with the r-process is to detect 
overabundances of Sr and Ba in the spectra of individual events as in the case of SN 1987A (see Sec. 
I4.3.3p . The simultaneous observations of light curves may also provide valuable information on the 
relationship between Fe production and r-process nucleosynthesis in core-collapse supernovae (see Sec. 

Em- 

The r-process involves a strong interplay between nuclear physics and astrophysics. For example, 
the conditions in the neutrino-driven wind associated with core-collapse supernovae are determined by 
the characteristics of neutrino emission and the mass and the radius of the protoneutron star, all of 
which are affected by the properties of hot and dense matter. The amount and the neutron-richness 
of the ejecta from neutron star mergers also depend on the behavior of matter at high density. Thus, 
the equation of state at high density plays a crucial role in both core-collapse supernova and neutron 
star merger models of the r-process. This equation of state also affects the gravitational radiation from 
neutron star mergers |167j . and therefore, may be constrained by the gravitational waves to be detected 
by experiments such as LIGO. At present, almost all of the essential nuclear physical input for the 
r-process reaction network has to be calculated from theory due to the lack of experimental data. This 
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situation may be changed by the Rare Isotope Accelerator (RIA), which is planned as the next major 
facility for nuclear physics in the United States |lti8j . Similar experimental efforts at RIKEN in Japan 
and GSI in Germany are also being discussed. Some of the issues regarding the r-process that may be 
addressed by RIA and similar facilities are considered below. 

The assumption of (n, 7) ^ (7, n) equilibrium greatly simplifies the reaction network for the r- 
process and provides a useful framework to discuss this process. It would be of great help to r-process 
calculations for different astrophysical environments if the conditions for (n, 7) ^ (7,"^) equilibrium 
(see Sec. 11.2. ip can be established by experimental data. The masses, the /3-decay lifetimes, and the 
neutron-capture cross sections for neutron-rich nuclei far from stability are needed for this purpose and 
may be obtained from experiments at RIA and similar facilities. 

Calculations with constant JH] or time-dependent |21^ neutron number density and temperature 
showed that an overstrong = 82 shell causes severe underproduction of the nuclei with A ~ 110-126. 
Quenching the strength of the N = 82 shell greatly alleviates the problem |lb9j . As the A^ = 82 shell 
must be sufficiently strong to produce the r-process abundance peak at A = 130, quenching should be 
most effective far below A = 130. Thus, the deficiency at the upper end of A ~ 110-126 may require 
other remedies such as neutrino-induced neutron emission from the progenitor nuclei in the abundance 
peak a.t A = 130 (see Sec. I3.1.1|) . The exact form of quenching may be identified by measurements of 
masses for the nuclei with A ~ 110-126 and A^ < 82 at RIA and similar facilities. 

Meteoritic data and the observed r-patterns in the ultra-metal-poor stars CS 22892-052 and CS 
31082-001 suggest that the high-frequency r-process events should produce the nuclei on both sides 
of A ~ 130 but very little of those with A ~ 130 (see Sec. E}. As a plausible scenario to explain 
this peculiar yield pattern, it was assumed that the r-process produces a freeze-out pattern covering 
190 < A < 320 with a peak at A ~ 195. Neutrino-induced fission of the progenitor nuclei during decay 
towards stability then produces the nuclei with 8Q < A < 190 but very little of those with A ~ 130 due 
to the characteristics of fission yields (see Sec. I4.3.2j) . The competition between neutron emission and 
fission in the deexcitation of unstable neutron-rich nuclei and the fission yields of such nuclei are crucial 
to the above scenario. Experiments at RIA and similar facilities may shed some light in this regard. 

A full understanding of the r-process requires quantitative description of the astrophysical envi- 
ronment, reliable nuclear physical input, and correct interpretation of meteoritic data and stellar ob- 
servations. If the role of neutrinos in the r-process can be established, then neutrino properties are 
an essential part of the problem. In view of the planning for measurements of nuclear properties at 
RIA and similar facilities, the ongoing observations of metal-poor stars with large telescopes (e.g., the 
Hubble Space Telescope, the Keck Telescopes, and the Very Large Telescope), and the ongoing and 
planned experiments on neutrino oscillations (e.g.. Super Kamiokande, Sudbury Neutrino Observatory, 
KamLAND, Mini-BooNE, and MINOS), much further progress in the understanding of the r-process 
is expected in the future. 
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